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The pseudogap Kondo problem, describing quantum impurities coupled to fermionic quasiparticles 
with a pseudogap density of states, p(a;) cx , shows a rich zero-temperature phase diagram, with 
different screened and free moment phases and associated transitions. We analyze both the particle- 
hole symmetric and asymmetric cases using renormalization group techniques. In the vicinity of 
r = 0, which plays the role of a lower-critical dimension, an expansion in the Kondo coupling 
is appropriate. In contrast, r = 1 is the upper-critical dimension in the absence of particle-hole 
symmetry, and here insight can be gained using an expansion in the hybridization strength of the 
Anderson model. As a by-product, we show that the particle-hole symmetric strong-coupling fixed 
point for r < 1 is described by a resonant level model, and corresponds to an intermediate-coupling 
fixed point in the renormalization group language. Interestingly, the value r — 1/2 plays the role of a 
second lower-critical dimension in the particle-hole symmetric case, and there we can make progress 
by a novel expansion performed around a resonant level model. The different expansions allow a 
complete description of all critical fixed points of the models and can be used to compute a variety of 
properties near criticality, describing universal local-moment fluctuations at these impurity quantum 
phase transitions. 



I. INTRODUCTION 

Non-trivial fixed points and associated phase transi- 
tions in quantum impurity problems have been subject 
of considerable interest in recent years, with applications 
for impurities in correlated bulk systems, in transport 
through nanostructures, and for strongly correlated lat- 
tice models in the framework of dynamical mean-field 
theory. Many of those impurity phase transitions occur 
in variations of the well-known Kondo model^i which de- 
scribes the screening of localized magnetic moments by 
metallic conduction electrons. A paradigmatic example 
of an intermediate-coupling impurity fixed point can be 
found in the two-channel Kondo effect. 

Non-metallic hosts, where the fermionic bath density 
of states (DOS) vanishes at the Fermi level, offer a differ- 
ent route to unconventional impurity physics. Of partic- 
ular interest is the Kondo effect in so-called pseudogap 
systemsSiSiliSiSiLS, where the fermionic bath density of 
states follows a power law at low energies, p{uj) oc iVo|a^|'' 
(r > 0). Such a behavior arises in semimetals, in cer- 
tain zero-gap semiconductors, and in systems with long- 
range order where the order parameter has nodes at the 
Fermi surface, e.g., p- and d-wave superconductors {r — 2 
and 1). Indeed, in d-wave high- Tc superconductors non- 
trivial Kondo-like behavior has been observed associated 
with the magnetic moments induced by non-magnetic Zn 
impurities^' Note that the limit r oo corresponds 
to a system with a hard gap. 

The pseudogap Kondo problem has attracted substan- 
tial attention during the last decade. A number of 
studiesSi^ employed a slave-boson large-A^ technique; 
significant progress and insight came from numerical 
renormalization group (NRG) calculations^iSii and the 
local moment approach^. It was found that a zero- 
temperature phase transition occurs at a critical Kondo 



coupling, Jc, below which the impurity spin is unscreened 
even at lowest temperatures. Also, the behavior depends 
sensitively on the presence or absence of particle-hole (p- 
h) asymmetry, which can arise, e.g., from a band asym- 
metry at high energies or a potential scattering term at 
the impurity site. A comprehensive discussion of possi- 
ble fixed points and their thermodynamic properties has 
been given in Ref. 6 based on the NRG approach. 

Until recently, analytical knowledge about the critical 
properties of the pseudogap Kondo transition was lim- 
ited. Previous works employed a weak-coupling renor- 
malization group (RG) method, based on an expansion 
in the dimensionless Kondo coupling j — iVo Jk. It was 
found that an unstable RG fixed point exists at j = r, 
corresponding to a continuous phase transition between 
the free and screened moment phaseA Thus, the per- 
turbative computation of critical properties within this 
approach is restricted to small r. Interestingly, the NRG 
studies^ showed that the fixed-point structure changes at 
r = r* « 0.375 and also at r=^, rendering the relevant 
case of 7' = 1 inaccessible from weak coupling. In the p-h 
symmetric case, for r > ^ the phase transition was found 
to disappear, and the impurity is always unscreened in- 
dependent of the value of Jk- In contrast, in the asym- 
metric case the phase transition is present for arbitrary 
r > 0. Numerical calculations^'^ indicated that the criti- 
cal fluctuations in the p-h asymmetric case change their 
character at r = 1: whereas for r < 1 the exponents take 
non-trivial r-dependent values and obey hyperscaling, ex- 
ponents are trivial for r > 1 and hyperscaling is violated. 
These flndings suggest to identify r = 1 as upper-critical 
"dimension" of the problem, whereas r — plays the role 
of a lower-critical "dimension". As known, e.g., from the 
critical theory of magnetsii, the description of the tran- 
sitions using perturbative RG requires different theoreti- 
cal formulations near the upper-critical and lower-critical 



2 



dimensions, i.e., the 0^ theory and the non-Unear sigma 
model in the magnetic case. 

In this paper, we provide a comprehensive analytical 
account of the phase transitions in the pseudogap An- 
derson and Kondo models, including the proper theories 
for the critical "dimensions". This is made possible by 
working with the Anderson instead of the Kondo model 
- the degrees of freedom of the Anderson model turn out 
to provide a natural description of the low-energy physics 
at the quantum phase transitions near r = ^ as well as 
close to and above r = 1. We shall consider epsilon- 
type expansions in the hybridization, the on-site energy, 
and the interaction strength. Those expansions lead to 
different theories for the p-h symmetric and asymmet- 
ric cases. Interestingly, in the pseudogap Kondo model 
the phase transitions near the lower-critical and upper- 
critical dimension are not adiabatically connected, as the 
fixed point structure changes both a,t r — r* and r — ^. 
Thus the present quantum impurity problem has a more 
complicated flow structure than the critical theory of 
magnets, where the (2 -I- e) and (4 — e) expansions are 
believed to describe the same critical fixed point. 

In the p-h symmetric case of the pseudogap Kondo 
problem the line of non-trivial phase transitions termi- 
nates at two lower- critical dimensions (!), r = and 
r = i. Near r = ^ we find an expansion around a non- 
interacting resonant level model, together with pertur- 
bative RG, to provide access to the critical fixed point, 
with the expansion being controlled in the small param- 
eter (i — r), see Sec. lIVI Interestingly, the weak-coupling 
expansion for the Kondo model, presented in Sec. IIIII 
provides a different means to access the same critical 
fixed point, but with the small parameter being r; the 
two expansions can be expected to match. 

In the p-h asymmetric case an expansion can be done 
in the hybridization around the valence-fluctuation point 
of the Anderson model. Bare perturbation theory is suffi- 
cient for all r > 1; for r < 1 a perturbative RG procedure 
is required to calculate critical properties, with the ex- 
pansion being controlled in the small parameter (1 — r). 
In particular, this identifies r = 1 as the upper-critical di- 
mension of the (asymmetric) pseudogap Kondo problem, 
and consequently observables acquire logarithmic correc- 
tions for r = 1. A brief account on the p-h asymmetric 
case and the expansion around r = 1 has been given in a 
recent paperi-. We note that the flow of the asymmetric 
Anderson model in the metallic case, r = 0, was discussed 
by Haldanei^: here all initial parameter sets with finite 
hybridizations flow towards the strong-coupling (singlet) 
fixed point. 

For all cases listed above, we show that the critical 
properties of the Anderson and Kondo models are identi- 
cal, and we calculate various observables in renormalized 
perturbation theory. To label the fixed points, we will 
follow the notation of Ref. 6. 

Before continuing, we emphasize that standard tools 
for metallic Kondo models, such as bosonization, Bethe 
ansatz, and conformal field theory, are not easily applica- 



ble in the present case of a pseudogap density of states, as 
the problem cannot be described using linearly dispersing 
fermions in one dimension. Furthermore, integrating out 
the fermions from the problem, in order to arrive at an 
effective statistical mechanics model containing impurity 
degrees of freedom only, cannot be performed easily: the 
fermionic determinants arising in (l-t-r) dimensions can- 
not be simply evaluated. This implies that the pseudogap 
Kondo model does not map onto a one-dimensional (e.g. 
Ising) model with long-ranged interactions, in contrast 
to e.g. the spin-boson modcl-'^''. Indeed, the phase transi- 
tions in the pseudogap Kondo model and the sub-ohmic 
spin-boson model are in different universality classesi^. 
Therefore we believe that our combined RG analysis pro- 
vides a unique tool for analyzing the pseudogap Kondo 
problem. 

A. Models 

The starting point of our discussion will be the single- 
impurity Anderson model with a pseudogap host density 
of states, H = Ha + H-h- 

Ha = eo/t./. + C/ori/T"/|+^o(/>.(0) + h.c.) , (1) 

.A 

Hh = dk \k\'' kcl^Cka 

J -A 

where we have represented the bath, Tib, by linearly 
dispersing chiral fermions Cka, summation over repeated 
spin indices a is implied, and Co-(O) = / dk\kYckcr is the 
conduction electron operator at the impurity site. The 
spectral density of the Cct(O) fermions follows the power 
law 1^1'' below the ultra-violet (UV) cutoff A; details of 
the density of states at high energies are irrelevant for 
the discussion in this paper. The four possible impurity 
states will be labelled with | T), I i) for the spin-carrying 
states, |c) for the empty and |d) for the doubly occupied 
state. Provided that the conduction band is p-h symmet- 
ric, the above model obeys p-h symmetry for Uq = —2eo 
- this p-h symmetry can be considered as SU(2) pseu- 
dospin, i.e., the full symmetry of the model is SU(2)spin 
X SU(2)charge- Asymmetry of the high-energy part of the 
conduction band has the same net effect as asymmetry of 
the impurity states; we will always assume that the low- 
energy part of the band is asymptotically symmetric, i.e., 
the prefactor of in the DOS is equal for positive and 
negative u. 

The transformation 

f(7 ^ fa 5 

Cka - (2) 

converts all particles into holes and vice versa, formally 
Eq —{sq -\- Uq), Vq —Vq. Physically, the roles of the 
states I t) and | |) are interchanged, as well as the states 
|e) and |d). It is useful to consider another transforma- 
tion, 

fl ^ fu /i — ^ /i ' 
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(3) 



which transforms | f) ^ |d), | l) ^ |e). Here, the spinful 
doublet of impurity states is transformed into the spinless 
doublet and vice versa, i.e., the two SU(2) sectors are 
interchanged. 

In the so-called Kondo limit of the Anderson model 
charge fluctuations are frozen out, and the impurity site 
is mainly singly occupied. Via Schrieffcr- Wolff transfor- 
mation one obtains the standard Kondo model, Ti. ~ 
Ti-K + Tih, with 



Hk = JkS ■ s(0) 



(4) 



where the impurity spin S is coupled to the conduction 
electron spin at site 0, Sq(0) — cj. (0)ct"^,Co-'(0)/2, and 
cr" is the vector of Pauli matrices. The Kondo coupling 
is related to the parameters of the Anderson model ^ 
through: 



Jk = 2 V 



1 

N 



1 



\Uo + eo\ 



(5) 



The Kondo limit is reached by taking Uq oo, Eq 
—oo, Vq oo, keeping Jk fixed. In the absence of p-h 
symmetry the Schrieffer- Wolff transformation also gen- 
erates a potential scattering term in the effective Kondo 
model^. 

In the absence of an external magnetic field all above 
models preserve SU(2) spin symmetry. Spin anisotropies 
turn out to be irrelevant at the critical fixed points, see 
Appendix^ The effect of a magnetic field will be briefly 
discussed in Sec. I Villi 
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B. Summary of results 

Our main results are summarized in the RG flow dia- 
grams in Figs. Hand 12 for the p-h symmetric and asym- 
metric cases, respectively. 

In the symmetric case, the ranges of exponent values 
r = 0, 0<r<i, i<r<l, and r > I lead to quite dif- 
ferent behavior, and are shown separately in Fig. ^ No 
transition occurs for r — 0: for any non-zero hybridiza- 
tion the flow is towards the metallic Kondo-screened fixed 
point (SC). This well-known fixed point can be identified 
as the stable fixed point of a resonant level model; we ar- 
gue below that this is actually an intermediate-coupling 
fixed point. 

For < r < i, small values of the hybridization leave 
the impurity spin unscreened provided that Eq < 0, i.e., 
there is a stable local-moment fixed point (LM) corre- 
sponding toe = — oo,w = 0. A transition line at negative 
e, with an unstable fixed point (symmetric critical, SCR) 
at finite v, |e|, separates the flow towards LM from the 
flow to the symmetric strong-coupling fixed point (SSC). 
The strong-coupling fixed point displays its intermediate- 
coupling properties now in a finite residual entropy and a 
finite magnetic moment, see Sec. lIVBl As r the SCR 



FIG. 1: Schematic RG flow diagrams for the particle-hole 
symmetric single-impurity Anderson model with a pseudogap 
DOS, p{i^) oc |tj|' . The horizontal axis denotes the renormal- 
ized on-site level energy e (related to the on-site repulsion u 
hy u = — 2e), the vertical axis is the renormalized hybridiza- 
tion V. The thick lines correspond to continuous boundary 
phase transitions; the full (open) circles are stable (unstable) 
fixed points, for details see text. All fixed points at non- 
zero e have a mirror image at —e, related by the particle-hole 
transformation @. a) r — Q, i.e., the familiar metallic case. 
For any finite v the flow is towards the strong-coupling fixed 
point (SC), describing Kondo screening, b) < r < i: The 
local-moment fixed point (LM) is stable, and the transition 
to symmetric strong coupling (SSC) is controlled by the SCR 
fixed point. For r — > 0, SCR approaches LM, and the critical 
behavior at SCR is accessible via an expansion in the Kondo 
coupling j. In contrast, for — > |, SCR approaches SSC, 
and the critical behavior can be accessed by expanding in the 
deviation from SCR, i.e., in e = —u/2. c) i < r < 1: n is 
still relevant at u = 0. However, SSC is now unstable w.r.t. 
finite u. At finite v, the transition between the two stable 
fixed points LM and LM' is controlled by SSC (which is now 
a critical fixed point!), d) r > 1: « is irrelevant, and the only 
transition is a level crossing (with perturbative corrections) 
occurring at u = u = 0, i.e., at the free-impurity fixed point 
(FImp). 
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fixed point merges with LM, in a manner characteristic 
for a lower-critical dimension, i.e., with diverging correla- 
tion length exponent. A second critical fixed point SCR' 
exists for e > which separates the symmetric strong- 
coupling phase (SSC) from one with a free charge doublet 
(LM'). 

As r ^ ^ the symmetric critical fixed points merge 
with the strong-coupling one, again in a manner charac- 
teristic for a lower-critical dimension. For r > ^ the fixed 
points SCR and SCR' cease to exist; the strong-coupling 
SSC fixed point becomes infrared unstable, and controls 
the LM - LM' transition. Finally, the structure of the 
flow changes again at r = 1: for r — > 1 the unstable 
strong-coupling flxed point (SSC) moves towards v = 0, 
i.e., the free- impurity fixed point (FImp), and for r > 1 
no non-trivial fixed point remains. 

For maximal p-h asymmetry, realized in the Anderson 
model through Uq — oo, one has to distinguish exponent 
ranges r = 0, < r < r* , r* < r < 1, and r > 1. In 
the metallic case r = any non-zero hybridization gen- 
erates fiow to strong coupling with complete screening 
- the strong-coupling fixed point is the same as in the 
p-h symmetric situation, as p-h symmetry is marginally 
irrelevant at strong coupling. For all r > the situation 
is drastically different: small Vq leaves the moment un- 
screened, whereas large Vq directs the flow towards a new, 
p-h asymmetric, strong-coupling fixed point (ASC). The 
character of the critical fixed point separating the two 
phases depends^ on r: for < r < r* p-h symmetry is 
restored, and the critical fixed point is the one of the p-h 
symmetric model. For r* < r < 1 there is a separate 
critical fixed point (ACR) which is p-h asymmetric, i.e., 
located at finite v and e. For r — > 1 the critical fixed point 
moves towards w — > 0, and for r > 1 the phase transition 
becomes a level crossing (with perturbative corrections) , 
controlled by the valence-fluctuation fixed point (VFl), 
see Fig. |21 

We finally discuss the general case of finite p-h asym- 
metry, more details will be given in Sec. I VII Power count- 
ing shows that LM (SSC) are always (un)stable w.r.t. p- 
h asymmetry. The symmetric critical SCR fixed point is 
stable w.r.t. p-h asymmetry for small r. In contrast, for 
r < i SCR is unstable towards p-h asymmetry, as it is 
close to SSC in this regime. This requires the existence of 
a specific r value where this change in character occurs: 
the is precisely r = r* 0.375 where p-h asymmetry 
at SCR is marginal^. Upon increasing r beyond r* the 
p-h asymmetric critical fixed point (ACR) splits off from 
SCR. In other words, upon approaching r* from large r 
the ACR fixed point moves towards small effective p-h 
asymmetry, and at r = r* ACR merges into SCR, im- 
plying p-h symmetry is dynamically restored. As stated 
above, the description of ACR using an expansion around 
VFl consequently breaks down as r ^ r*"*". Neither from 
numerics^ nor from the present RG are there indications 
for the existence of a second asymmetric critical fixed 
point besides ACR; thus, the critical properties for finite 
p-h asymmetry are always equivalent to the ones of a 



a) r* < r < 1 ^ A 




e 



FIG. 2: Schematic RG flow diagram for the maxi- 
mally particle-hole asymmetric pseudogap Anderson impu- 
rity model. The horizontal axis denotes the on-site impurity 
energy, e, the vertical axis is the fermionic coupling v, the 
bare on-site repulsion is fixed at wo = tx). The symbols are 
as in Fig. a) r* < r < 1: w is relevant, and the tran- 
sition is controlled by an interacting fixed point (ACR). As 
r ^ r* ~ 0.375, p-h symmetry at the critical fixed point is 
dynamically restored, and ACR merges into the SCR fixed 
point of Fig. Q - this cannot be described using the RG of 
Sec. [3 In the metallic r = situation, studied by Haldana— , 
the flow from any point with n 7^ is towards the screened 
singlet fixed point with £ = 00. b)r>l: wis irrelevant, and 
the transition is a level crossing with perturbative corrections, 
occuring at d = £ = 0, i.e., the valence-fluctuation fixed point 
(VFl). 

model with maximal p-h asymmetry. 

Taken together, the above observations show that r = 
plays the role of a lower-critical dimension: as r — O"*" , 
the correlation length exponent diverges, and the second- 
order transition turns into a Kosterlitz-Thouless transi- 
tion at r = 0. Interestingly, in the symmetric case the 
correlation length exponent also diverges as r ^ i , 
and the transition between LM and SSC disappears for 
r>i: r — ^ is a second lower-critical dimension for the 
p-h symmetric problem. In the asymmetric case, there 
is a transition between LM and ASC for all r > 0, and 
r = 1 is equivalent to the upper-critical dimension, above 
which the critical fixed point is non-interacting (actually 
a level crossing). 



C. Outline 

The rest of this paper is organized as follows: Sec. IH] 
introduces the observables to be evaluated in the course 
of the paper, together with their expected scaling be- 
havior near criticality. In Sec. IIIII we briefiy review the 
standard weak-coupling perturbative RG for the Kondo 
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model, which is suitable to describe the quantum phase 
transition for small r. Sec. |W] discusses the particle- 
hole symmetric Anderson model. Starting from the non- 
interacting case, £o = Uq = 0, we first discuss the 
physics of the resulting non-interacting resonant level 
model - interestingly this can be identified with a sta- 
ble intermediate-coupling fixed point. We then use a 
perturbative expansion in Uq to access the critical fixed 
points for r i. In Sec. we turn to the situation 
with maximal p-h asymmetry, i.e., Uq = oo, and show 
that an expansion in the hybridization provides access 
to the critical properties for r > 1 as well as for r ^ 1. 
In Sec. IVII we consider the case of general p-h asymme- 
try. Sec. IVIIII briefly describes the effect of a magnetic 
field: the pseudogap model is shown to permit a sharp 
transition as function of a field applied to the impurity 
for couplings larger than the zero-field critical coupling. 
In Sec. IVIII wc compare the physics of the Anderson and 
Kondo models, arguing that the transitions in both mod- 
els fall in the same universality classes. A brief discus- 
sion of applications concludes the paper. All renormal- 
ization group calculations will employ the field-theoretic 
RG scheme together with dimensional regularization 
and minimal subtraction of poles, with details given in 
the appendices; one-loop RG results can equivalently be 
obtained using the familiar momentum-shell method. 

II. OBSERVABLES AND SCALING 

To establish notations and to pave the way for the RG 
analysis below, we introduce a few observables together 
with their expected scaling properties. 

A. Susceptibilities 

Magnetic susceptibilities are obtained by coupling an 
external magnetic field to the bulk electronic degrees of 
freedom in Tib, 

-iJua(2:)(4^".'C.')(^) (6) 
and to the impurity part TYa, ^k, 

Himp.aif^CTfyfy' fa') j -ffimp.Q'S'a (7) 

for the Anderson Q and Kondo models, respectively. 
The bulk field varies slowly as function of the space 
coordinate, and i/i,np is the magnetic field at the location 
of the impurity. 

With these definitions, a spatially uniform field applied 
to the whole system corresponds to ~ i?imp = H. 
Response functions can be defined from second deriva- 
tives of the thermodynamic potential, = —TlnZ, in 
the standard waj*ii: Xu,u measures the bulk response to 
a field applied to the bulk, Ximp,imp is the impurity re- 
sponse to a field applied to the impurity, and Xu,imp is 
the cross-response of the bulk to an impurity field. 



The impurity contribution to the total susceptibility is 
defined as 

Ximp,imp ~l" ijnp H" (Xu.u Xu,u (8) 

where Xu"u'^ is the susceptibility of the bulk system in ab- 
sence of the impurity. For an unscreened impurity spin 
of size 5 = i we expect Ximp(T ^ 0) = l/(4r) in the 
low-temperature limit, and this is precisely the result in 
the whole LM phase. A fully screened moment will be 
characterized by Tximp = 0; note that the SSC fixed 
point displays a finite value of Tximp for r > 0. At crit- 
icality Ximp does not acquire an anomalous dimension — 
(in contrast to xioc below), because it is a response func- 
tion associated to the conserved quantity ^tot- Thus we 
expect a Curie law 

^imximp(T) = ^, (9) 

where the prefactor Cimp is in general a non-trivial uni- 
versal constant different from the free-impurity value 
S{S + l)/3. Apparently, Eq. JSJ can be interpreted as 
the Curie response of a fractional effective spini^. 
The local impurity susceptibility is given by 

Xloc 

(T) 

Ximp, imp 7 

(10) 

which is equivalent to the zero-frequency impurity spin 
autocorrelation function. In the unscreened phase we 
have Xloc oc 1/T as T — > 0; we can consider this as arising 
from the overlap of the local impurity moment with the 
total, freely fluctuating, moment of 5 = 1/2, and so write 

2 

liinxioc(r) = ^. (11) 

The quantity mioc turns out to be a suitable order 
parameteriSii for the phase transitions between an un- 
screened and a screened spin: it vanishes continuously as 
t 0~ , where t is the dimensionless measure of the dis- 
tance to criticality; in the Kondo model t = ( Jk — Jc) / Jc, 
whereas in the Anderson model t = {Vo — Vqc)/Vo^. Thus, 
^Xioc is not pinned to the value of 1 /4 for t < (in con- 
trast to Tximp)- Remarkably, mjoc = at the SSC fixed 
point for r < 1, although Tximp = r/8 there. 

The phase transitions occurring for < r < 1 are de- 
scribed by interacting fixed points, and thus obey strong 
hyperscaling properties, including lo/T scaling in dynam- 
ical quantitiesSSi. For instance, the local dynamic suscep- 
tibility will follow a scaling form 

which describes critical local-moment fluctuations, and 
the local static susceptibility obeys 

Xioc(T) = ^<i>2(^^j . (13) 
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Here, rj^ is a universal anomalous exponent, which con- 
trols the anomalous decay of the two-point correlations of 
the impurity spin, and i>i_2 are universal crossover func- 
tions (for the specific critical fixed point and for fixed r) , 
whereas Bi^2 are non-universal prefactors. Furthermore, 
v is the correlation length exponent, describing the flow 
away from criticality: when the system is tuned through 
the transition, the characteristic energy scale T*, above 
which critical behavior is observed, vanishes as^S 



T* cx lir; 



(14) 



the dynamical critical exponent z can be set to unity in 
the present (0 + l)-dimensional problem. Note that at 
criticality, t = 0, the relation (|13|) reduces to xioc(T) « 

Hyperscaling can be used to derive relations between 
critical exponents. The susceptibility exponent ry^ and 
the correlation length exponent of a specific transi- 
tion are sufficient to determine all critical exponents as- 
sociated with a local magnetic field^. In particular, the 
T ^ local susceptibility away from criticality obeys 

X\oc{t > 0) cx t"^ , 7 = 1/ (1 - , 
Txioc(t <0) cc i-tr' , 7' = i^%, (15) 

which can be derived from a scaling ansatz for the im- 
purity part of the free energy^. The last relation implies 
the order parameter vanishing as 



mioc oc (-t)'^''-/2 , 



(16) 



Note that hyperscaling holds for all critical fixed points 
of the pseudogap Kondo problem with < r < 1. 



B. Impurity entropy 

In general, zero-temperature critical points in quan- 
tum impurity models can show a finite residual entropy 
[in contrast to bulk quantum critical points where the 
entropy usually vanishes with a power law, S{T) cx T^]. 
For the models at hand, the impurity contribution to the 
low-temperature entropy is obtained by a perturbative 
evaluation of the thermodynamic potential and taking 
the temperature derivative. This will yield epsilon-type 
expansions for the ground state entropy S'imp(T — 0), 
with explicit results given below. 

Note that the impurity part of the thermodynamic po- 
tential will usually diverge with the cutoff, i.e., we have 
^^imp = ^'imp " ^''S'imp, whcrc E'imp is the non- universal 
(cutoff-dependent) impurity contribution to the ground- 
state energy. However, the impurity entropy 5i,np is fully 
universal, and the UV cutoff can be sent to infinity after 
taking the temperature derivative of llimp- 

Thermodynamic stability requires that the total en- 
tropy of a system decreases upon decreasing tempera- 
ture, dTS{T) > 0. This raises the question of whether 
the impurity part of the entropy. Simp has to decrease 



under RG flow (which is equivalent to decreasing T). 
The so-called g-theorem.^^ provides a proof of this conjec- 
ture for systems with short-ranged interactions; for most 
quantum impurity problems this appears to apply. In- 
terestingly, the pseudogap Kondo problem provides an 
explicit counter-example, as the two critical fixed points 
obey S'scR < 5'acr, with the RG flow being from SCR 
to ACR (!), see Sec.0for details. (For another counter- 
example see Ref. li^) 



C. T matrix 

An important quantity in an Anderson or Kondo model 
is the conduction electron T matrix, describing the scat- 
tering of the c electrons off the impurity. For an Anderson 
model, the T matrix is just given by T{uj) = VQGf{uj) 
where G/ is the full impurity / electron Green's function. 
For a Kondo model, it is useful to deflne a propagator 
Gt of the composite operator = f^fa'Ca-', such that 
the T matrix is given by T(w) = J^Gt(<^)- As with the 
local susceptibility, we expect a scaling form of the T ma- 
trix spectral density near the intermediate-coupling fixed 
points similar to Eq. (|12f) . In particular, at criticality a 
power law occurs: 



r(c.) 



(17) 



Remarkably, we will find the exact result rjT — 1 — r for 
r < 1, i.e., for all interacting critical fixed points con- 
sidered in this paper the T matrix follows T{uj) cx lu^^. 
NRG calculations have found precisely this critical diver- 
gence, for both the symmetric and asymmetric critical 
points^. At the trivial fixed points (LM, ASC) the be- 
havior of the T matrix follows from perturbation theory 
in the hybridization to be T{uj) cx lu^ . 

Notably, the T matrix can be directly observed 
in experiments, due to recent advances in low- 
temperature scanning tunneling microscopy, as has 
been demonstrated, e.g., with high-temperature 
superconductor322iSii2^. 



D. Phase shifts 

Fixed points which can be described in terms of free 
fermions can be characterized by the s-wave conduc- 
tion electron phase shift, 6q{lu), which can be related 
to the conduction electron T matrix through ^o('-^) — 
argT(cj). A decoupled impurity simply has a phase shift 
do — 0, whereas a p-h symmetric Kondo-screened impu- 
rity in a metallic host has a low-energy phase shift of 
So{uj) — ■|sgn(— tj). A detailed discussion for the pseu- 
dogap model has been given in Ref. 6, in the body of the 
paper we will simply quote the results. 
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III. WEAK-COUPLING RG FOR THE KONDO 
MODEL 

In this section we briefly summarize the weak-couphng 
RG for the pseudogap Kondo model Q), as first discussed 
by Withoff and Fradkin^. Perturbative RG is performed 
around Jk = 0, i.e., the local-moment fixed point (LM): 
this will allow to access the (p-h symmetric) critical fixed 
point SCR, which is located close to LM for small DOS 
exponents r. 

A. Lower critical dimension: Expansion around the 
local-moment fixed point 

The RG flow equation for the renormalized Kondo cou- 
pling j, to two-loop order, read^ 

/3(j) = n-/ + y- (18) 
This yields an infi'ared unstable fixed point at 

J*=r+y+0(r3) (19) 

which controls the transition between the decoupled LM 
and the Kondo-screened SSC phases. The small-j expan- 
sion (|18|l - which is nothing but the generalization of An- 
derson poor man's scaling^'' to the pseudogap case - can- 
not give information about the strong-coupling behavior, 
and it can only describe critical properties for small r. 
(In the p-h symmetric case, the fixed point structure does 
not change within the exponent range < r < i, thus 
the present expansion is in principle valid up to r = i.) 

Adding a potential scattering term Vb gives a finite p-h 
asymmetry. Under RG, we find that V renormalizes to 
zero for r > 0, I3{V) = rV. Thus, within the range of 
applicability of the weak-coupling RG, p-h asymmetry is 
irrelevant. (Strictly, this applies for r < r* , see Ref. 6.) 

B. Observables near criticality 

We quote a few properties of the critical regime which 
have been determined in Ref. l26i Expanding the beta 
function 1)18(1 around the fixed point value H19() yields the 
correlation length exponent v. 

i=r-y + 0(r3). (20) 

The low-temperature impurity susceptibility and entropy 
at criticality are given by 

TXirap = i(l-r)+0(r2), (21) 

5i„,p = \n2(l + ^ A +0{r'). (22) 



The anomalous exponent of the local susceptibility eval- 
uates to 

Vx = r' + 0{r') . (23) 

A comparison of the above results with numerical data 
is given in Figs.OlElEl and[3below. 

Most importantly, the continuous transition controlled 
by the fixed point (|19|l . which exists only for r > 0, 
evolves smoothly into the Kosterlitz-Thouless transition 
at r = 0, J = 0, which separates the antiferromagnetic 
and the ferromagnetic metallic Kondo model. This is also 
indicated by the divergence of the correlation length ex- 
ponent (|20|l as r ^ 0+. Thus, r = can be identified as a 
lower-critical "dimension" of the pseudogap Kondo prob- 
lem. It is interesting to compare the present expansion 
with the (2-|-e) expansion for the non-linear sigma model, 
appropriate for magnets close to the lower-critical dimen- 
sion. The expansion is done about the ordered magnet, 
thus the LM phase with In 2 residual entropy takes the 
role of the ordered state in the pseudogap Kondo prob- 
lem. 



IV. PARTICLE-HOLE SYMMETRIC 
ANDERSON MODEL 

In the following sections of the paper we shift our at- 
tention from the Kondo model to the impurity Anderson 
model with a pseudogap density of states. This formula- 
tion will provide new insights into the RG fiow and the 
critical behavior of both the Anderson and Kondo mod- 
els. 

The coupling between impurity and bath is now the 
Anderson hybridization term, Vq, which turns out to be 
marginal in a RG expansion around Vq — for the DOS 
exponent r = 1 (in contrast to the Kondo coupling Jk 
which is marginal for r = 0). As we will show in Sec. Ivl 
the Anderson model provides the relevant low-energy de- 
grees of freedom for the p-h asymmetric pseudogap tran- 
sition near its upper-critical dimension. 

Interestingly, also the p-h symmetric version of the 
Anderson model allows to uncover highly non-trivial 
physics, in particular the special role played by the DOS 
exponent r = ^, where the transition disappears in the 
presence of p-h symmetry. Thus we start our analysis 
with the particle- hole symmetric Anderson model l(T]l. 
i.e., we keep Uq — — 2eo and discuss the physics as func- 
tion of Vo and eo- 



A. Trivial fixed points 

For vanishing hybridization, Vq = 0, the symmetric 
Anderson model |^ features three trivial fixed points: 
for £o < the ground state is a spinful doublet - this 
represents the local- moment fixed point (LM). For eo > 
we find a doublet of states (empty and doubly occupied). 
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denoted as LM' and related to LM by the p-h transfor- 
mation Both LM and LM' have a residual entropy of 
•Simp = ln2. At £o = a level crossing between the two 
doublets occurs, i.e., all four impurity states are degen- 
erate - this is the free- impurity fixed point (FImp), with 
residual entropy In 4. The impurity spin susceptibilities 
are 



Tx 



imp 



1/4 LM 
1/8 FImp 
LM' 



(24) 



the conduction electron phase shift is zero at all these 
fixed points. The hybridization term, Vb, is irrelevant at 
both the LM and LM' fixed points for r > 0, whereas for 
r = it is marginally relevant, as shown by the RG in 

sec inra 



B. Resonant level model: Intermediate-coupling 
fixed point 

It proves useful to discuss the Eq = = Q case, i.e., 
the physics on the vertical axis of the flow diagrams in 
Fig. ^ This non-interacting system is known as resonant 
level model, as the two spin species are decoupled. The 
problem can be solved exactly: the / electron self-energy 



(25) 



where Gco is the bare conduction electron Green's func- 
tion at the impurity location R = 0. In the low-energy 
limit the / electron propagator is then given by 

Gf{iuj.ny^ = iuJn - iA(3 sgn(cj„) (26) 

where the \u}nY self-energy term dominates for r < 1, 
and the prefactor Aq is 



(27) 



Before stating results for observables it is interesting 
to tackle the problem using RG techniques, with an ex- 
pansion in the hybridization strength Vq around the free- 
impurity fixed point (FImp, Vq = 0). We study the action 



Jo 



dr 



[fadrU+Vo{f.C,{0)+C 
A 

dk\k\'^ Ckaidr - k)cka 



(28) 



where (0) is the bath fermion field at the impurity po- 
sition as above. Power counting w.r.t. the Vb = fixed 
point, using dim[/] = 0, dim[c(0)] = (1 + r)/2, yields 



dim[Vb] 



l-r 



(29) 



i.e., the hybridization is relevant only for r < 1. 

To perform RG within the field-theoretic scheme^Si, we 
introduce a renormalized hybridisation v according to 
Vb = (Z^[f I \fZi)v^ where ^ is a renormalization energy 
scale, and and Z are the interaction and field renor- 
malization factors. The RG flow equation for v is found 
to be 



/3(w) = —TV -|- 



(30) 



Remarkably, this result is exact to all orders in per- 
turbation theory: the cubic term arises from the only 
self-energy diagram of the / fermions, and no vertex 
renormalizations occur (Zy = 1). This implies that the 
low-energy physics of the non-interacting resonant level 
model is controlled by the stable intermediate-coupling 
fixed point located at 



1 



(31) 



for < r < 1, which also applies to the familiar metallic 
case r = 0. 

The intermediate-coupling nature of the stable fixed 
point, with associated universal properties, is consistent 
with the results known from the exact solution of the 
problem, e.g., a universal conduction electron phase shift, 
a universal crossover in the temperature-dependent sus- 
ceptibility etc. 

We proceed with calculating a number of observables 
for the pseudogap resonant level model. Interestingly, 
this can be done in two ways: either (i) via the ex- 
act solution of the problem, i.e., by integrating out the 
c fermions exactly [leading to the propagator (|26|l ]. or 
equivalently (ii) by evaluating pcrturbative corrections to 
the FImp fixed point using the RG result (|31|l . utilizing 
standard renormalized perturbation theory, and noting 
that all corrections beyond second order in v vanish ex- 
actly within this scheme. Details of the calculation are 
in Appendix IXI 

We start with evaluating spin susceptibilities - note 
that we have kept two spin species in the model. The 
zero-temperature dynamic local susceptibility is propor- 
tional to the bubble formed with two / propagators H26|l . 



,l-2r 



(0 < r < 1) 



S{u;)^ (r>l) 



(32) 



where the case of r = i receives logarithmic corrections, 
see below. The low-temperature limit of the impurity 
susceptibility is found to be 



TXimpiT) 



the impurity entropy is 



Simp 2r In 2 , 



(33) 



(34) 



where the two last equations are valid for < r < 1; for 
r > 1 the resonant level model flows to the free-impurity 
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fixed point (FImp) with properties listed in Sec. IIV Al 
The conduction electron phase shift near the Fermi level, 
determined in Ref . 0, is 



sgn(-w) 



(1 -r)f (0 < r < 1) 

2 In I I ('^ = 1) 

0{uj'-^) (r > 1) 



(35) 



Interestingly, the resonant level model describes a 
screened impurity only in the metallic case. For the pseu- 
dogap case, r > 0, Eqs. (p?^ and show that the impu- 
rity is only partially screened: in a model of free fermions 
we have a residual entropy! We will see below that the 
resonant level model fixed point (|31|l can be identified 
with the symmetric strong-coupling fixed point (SSC) of 
Gonzalez-Buxton and Ingersent^ , introduced for the p-h 
symmetric Kondo and Anderson models. 



C. Expansion around the resonant level fixed point 

After having analyzed the behavior of the Anderson 
model in the non-interacting case, we proceed to study 
the stability of the resonant level fixed point w.r.t. a fi- 
nite interaction strength Uq, keeping p-h symmetry. Im- 
portantly, this fixed point, characterized by a finite hy- 
bridization strength between impurity and bath, is stable 
for small r [see Eq. (|37ll ]: we conclude that it can be iden- 
tified with the SSC fixed point of Ref. 6. 

Numerical resultsSife indicate that the quantum phase 
transition between LM and SSC disappears as r is in- 
creased to ^ , where the p-h symmetric critical fixed point 
(SCR) merges with the SSC fixed point. We shall show 
that an expansion around the SSC fixed point captures 
the physics of the SCR fixed point for f < i. Thus, this 
expansion describes the same critical fixed point as the 
weak-coupling expansion of Sec. IIIII but approaching it 
from r = ^ instead of r = 0. (The r values and i 
are two lower-critical dimensions for the p-h symmetric 
pseudogap Kondo problem.) 

The RG expansion below will be performed around 
an intermediate- coupling fixed point, in contrast to most 
analytical RG calculations which expand around trivial 
(i.e. weak or strong-coupling) fixed points. Strategically, 
one could think about a double expansion in Vq and Uq. 
However, this is not feasible, as the marginal dimensions 
for both couplings are different, r — 1 and r = ^, respec- 
tively. Therefore we choose to first integrate out the c 
fermions exactly, and then use standard RG tools for the 
expansion in Uo- 

Consequently, the starting point is the action 



S = ^ /^(w„) [iAoSgn(tj„)|tj„|''] /a.(w„) (36) 

where the / fermions are now "dressed" by the conduc- 
tion lines. ^0 is the non-universal number given in Eq. 



(|27|l . and we have assumed < r < 1. The interac- 
tion term has been written in a p-h symmetric form. A 
note is in order regarding the cutoff: The original model 
had a UV cutoff A, and this sets the upper bound for the 
\uJn\^ behavior of the / propagator (|26() . i.e., A is now the 
energy cutoff for the spectral density of the / fermions. 
The RG to be performed below can be understood as 
progressive reduction of this cutoff (although we will use 
the field-theoretic scheme where the cutoff is implicitely 
sent to infinity at an early stage). 

Dimensional analysis w.r.t. the Uq — situation, using 
dim[/] = (1 — r)/2, results in 



dim[[/o] = 2r - 1 = -e , 



(37) 



hence, for r > 1/2 the interaction term is relevant and 
the SSC fixed point is unstable. 

We continue with the RG analysis of (|36|l . To perform 
a perturbative expansion in Uq using the field-theoretic 
schemei^, we introduce a renormalized field and a dimen- 
sionless coupling according to 



fa = VZfR, 
12 



Un = 



Z2 



(38) 
(39) 



where /i is a renormalization energy scale as usual, and 
(— e) is the bare scaling dimension of Uq; we have ab- 
sorbed the non-universal number appearing in the 
dynamic term of the action (|36|l in order to obtain a uni- 
versal fixed point value of u. 

The complete RG analysis, needed to determine the 
flow of u, is presented in Appendix ^ here we restrict 
ourselves to the final results. P-h symmetry prohibits 
the occurence of even powers of u in the beta function 
of u, therefore the lowest contributions arise at two- loop 
order. Remarkably, no singular propagator renormaliza- 
tions occur, thus 



Z = 1 



(40) 



to all orders in perturbation theory. The RG flow equa- 
tion for the renormalized interaction u, arising now only 
from two-loop vertex renormalizations, is found to be 



3(7r - 2 In 4) 



(41) 



Next-to-leading order contributions would require a four- 
loop calculation which we do not attempt here. For pos- 
itive e, i.e., r < i, Eq. yields a pair of unstable 
fixed points at finite |u| (in addition to the stable one at 
u = 0); the correlation length exponent of the transition, 

Eq. H43I) . diverges as r ^ ^ . Thus, the behavior below 
r = i is similar to the standard behavior above a lower- 
critical dimension (e.g. in the non-linear sigma model for 
bulk magnet case). 
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D. 



r = 



Clearly, the unstable finite-u fixed points predicted by 
the perturbative RG equation H41|) do not necessarily ex- 
ist for the metallic case r = 0, as e = 1 is possibly out- 
side the convergence radius of the expansion. Indeed, the 
Kondo RG of Sec. 1111 Al shows that, at r = 0, the LM and 
LM' fixed points are unstable w.r.t. finite impurity cou- 
pling. As the resonant-level fixed point at u = is sta- 
ble, we conclude that the flow is directly from LM (LM') 
to SC, and SC represents the familiar strong-coupling 
Kondo fixed point, with complete screening of the spin. 
The RG flow is in Fig. [T]i. 



E. 



< r < 







For r values smaller than i, both the v — v" ,u 
flxed point (SSC) and the v — 0, \e\ = oo fixed points 
(LM, LM') are stable, and should be separated by criti- 
cal fixed points. The RG equation H41(l yields a pair of 
infrared unstable fixed points at 



3(7r-21n4) 



(42) 



with e = 1 — 2r. These two fixed points represent SCR 
and SCR', see the fiow diagram in Fig.^. Note that p-h 
symmetry also dictates that the fiow trajectories out of 
the SSC fixed point are horizontal in the u-v^ diagram. 
Therefore, close to r = i the SCR and SCR' fixed points 
are completely described by the fixed point coupling val- 
ues V* ^ and u* 



< r < 1 



1 r 

0.8 : 
0.6 : 

0.4 - 
0.2 - 







o ° ° □ V 



0.2 0.4 0.6 O.i 
r 



FIG. 3: Inverse correlation length exponent l/i/ obtained 
from NRG, at both the symmetric (squares) and asymmetric 
(triangles) critical points, together with the analytical RG 
results from the expansions in r fSec. IIIll Eq. 12UII . solid], in 
(I - r) [Sec. |TV| Eq. dashed], and in (1 - r) [Sec. El 

Eq. 1721 . dash-dot]. The numerical data have been partially 
extracted from Ref. using hyperscaling relations; for the 
symmetric model data are from Ref. |5j- 



between the two stable fixed points LM and LM'. The 
resulting flow diagram is in Fig. 



H. Observables near criticality 

Here we discuss critical properties of the SCR fixed 
point, the properties of SCR' are identical when trans- 
lated from spin to charge degrees of freedom. The corre- 
lation length exponent follows from expanding the beta 
function H41() around its fixed-point value: 



1 



(43) 



For r > i (r = i) the self-interaction u is a A comparison with results from NRG is shown in Fig. 13 



(marginally) relevant perturbation at the SSC fixed 
point, and 141|l does not yield additional non-trivial fixed 
points. As LM and LM' are stable, we can conclude that 
the fiow is from SSC directly towards LM (LM') for pos- 
itive (negative) Uq. Hence, SSC has become a critical 
fixed point, controlling the transition between LM and 
LM', which occurs at Uq = for any finite Vq- We shall 
not consider this transition in greater detail, apart from 
stating its correlation length exponent, l/iy = —e. Fig.^ 
displays the flow diagram arising from this discussion, be- 
ing consistent with the numerical results of Ref. 0. 



Close to r = i, the analytical expression nicely matches 
the numerical results, however, higher-order corrections 
in the expansion quickly become important. 

We continue with the quantities introduced in Sec. 
As usual for an expansion where the non-linear coupling 
has an infrared unstable fixed point (as occurs above the 
lower-critical dimension in standard situations), the UV 
cutoff needs to be kept explicitly, and intermediate quan- 
tities will diverge with the UV cutoff (see e.g. the calcula- 
tion of the impurity entropy). However, these divergences 
will cancel in the final expressions for universal observ- 
ables, and this is an important check for the consistency 
of our calculations. 



G. 



r > 1 



The physics of the symmetric Anderson model for 
r > 1 is easily discussed: the hybridization term Vq is 
irrelevant for all Uq, Eq. (|29|l . The free-impurity fixed 
point (FImp) is the only remaining fixed point at Uq — 0. 
It is unstable w.r.t. finite Uq, and controls the transition 



1. Local susceptibility 

The local susceptibility at the SSC fixed point, i.e., at 
tree level, follows the power law xioc oc uj^~''^^ = w*^. To 
obtain corrections to the tree-level result, one introduces 
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a) 



b) 



c) d) 



e) 




FIG. 4: Feynman diagrams for the symmetric Anderson im- 
purity model, where the expansion is done in Uo around the 
resonant level model fixed point. Full lines are dressed fa 
propagators with the self-energy arising from the conduction 
electrons already taken into account, i.e., with the propagator 
given in Eq. 12611 . a) Bare interaction vertex Uo- b) Corre- 
lation function involving a composite (//) operator (wiggly 
line), together with the first perturbative correction, c) Uo 
correction to the local susceptibility Xi mp,inip — Xioc, opeu Cir- 
cles are sources, d) Uo correction to Xu,m which contributes 
to the Curie term of the impurity susceptibility Ximp. Dashed 
lines are conduction electron lines, and the full dots are the 
Vb vertices, e) Uq contribution to the impurity free energy. 



1 r- 

0.8 : 
0.6 : 

0.4 - 

0.2 - 

0^ 
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FIG. 5: NRG datai for the local susceptibility exponent, 
»7x, defined through xioc oc T~^+''x^ at both the symmetric 
(squares) and asymmetric (triangles) critical points, together 
with the renormalized perturbation theory results from the 
expansions in r [Sec. lIIII Eq. 12311 . solid], in ( | — r) [Sec. lIVI 
Eq. (EZl, dashed], and in (1 - r) [Sec.|V| Eq. JTU, dash-dot]. 



2. Impurity susceptibility 



a xioc renormalization factor, Z^, from which one obtains 
the anomalous exponent according to 



du 



(44) 



Note that 1 — % = e+^Z^' ^^^^ defined in Eqs. (|12I13|I . 
due to the non-trivial structure of the problem already 
at tree level. 

Different ways can be used to determine Z^. Realizing 
that xioc is a correlation function of a composite (//) 
operator leads to Z^ — Z\^ where Z2 is the renormaliza- 
tion factor associated with (//). Z^ can be calculated 
from the correlation function shown in Fig. which 
receives a perturbative correction to first order in C/q. 
Alternatively, xioc can be calculated directly, and a sin- 
gle diagram gives a contribution of order Uq (Fig. ^). 
Both ways lead to 



2u 

= 1 + — ■ 



(45) 



with details given in Appendix^ The result for xioc is 
proportional to the non- universal number A^"^ ^ however, 
the exponent is universal: 



7y' =-u* = 2[3(^-21n4)]-i/2^ 



(46) 



The local susceptibihty thus follows xioc(T) oc T ^^"^^ 
with 



r7x = 2 - 2r + 2.688^/^ - r 



(47) 



where the first term contains the tree level expression, 
and further 0(e) terms arise from higher-order pertur- 
bative corrections. A comparison with NRG results is 
given in Fig. |S1 where good agreement near r 
observed. 



i can be 



For the impurity contribution to the uniform suscep- 
tibility we expect a Curie law, as discussed in Sec. ^1 
At tree level, a term of Curie form does only arise from 
Xu,u, with Txu.u = r/8. Both Ximp,imp and Xu,imp are 
less singular for r < 1, consistent with Txioc = 0, see 
also Appendix Also note that Xu,u is the only of the 
three terms where the non-universal number ^0 drops 
out. 

We are interested in corrections to Tximp to lowest 
non-trivial order in Uq, and consequently those correc- 
tions can only occur in Xu,u- A single diagram contributes 
to first order in Uq (Fig. 0Ji) , which gives 



^0 



AXimp = f^O^ 



dklkl'-Tj^ 



{iujn - ky 



(48) 



The k integral can be performed first, with the UV cutoff 
sent to infinity, and the frequency summation then leads 
to 



Aximp = ^ (1 - 2" 



27r2(i+'-) 



^0 



(49) 



In the low-energy limit, the combination {UqT^/Aq) 
approaches a universal value, the non-universal prefac- 
tors Vq are seen to cancel, and the result has the expected 
universal Curie form. Introducing the renormalized cou- 
pling u we have to leading order in e: 



A(Txh„p) 



1 



1 



2V2 



C(3/2)^ 
27r3 



(50) 



Using the fixed point value of u (I42|) . the result for the 
impurity susceptibility at the SCR fixed point reads 



0.1942 



r + Oie) (51) 



12 



0.25 




Co 




FIG. 6: Numerical data for the impurity susceptibility, 
TXimp, at both the symmetric (squares) and asymmetric (tri- 
angles) critical points, together with the renormalized per- 
turbation theory results from the expansions in r [Sec. 11111 
Eq. (HIJ, solid], in (| - r) [Sec.|TVl Eq. (EU, dashed], and in 
(1 - r) [Sec. El Eq. (gKJ, dash-dot]. The thin line is the value 
of the SSC fixed point, Tyiinp = r/8. The numerical data are 
partially taken from Ref. |a; we have re-calculated the data 
points near r = 1, as the logarithmically slow flow at r = 1 
complicates the data analysis. 



FIG. 7: As Fig. |S| but for the impurity entropy Simp; the 
NRG values at the SSC fixed point are hardly distinguishable 
from In 2. The perturbative expansions are in Eq. 12211 (solid), 
Eq. I55II (dashed), and Eq. II89II (dash-dot); the thin line is 
the value of the SSC fixed point, Tximp = 2rln2. 



and T, with an extrapolation of the result to A/T — > oo, 
and obtained [9t-.] = 0.1590 ± 0.0005, i.e., 



where r/8 is the tree level contribution, and the 0{e) 
represents higher perturbative terms as above. This re- 
sult can be nicely compared to NRG data of Ref. |6|, see 
Fig.El 



3. Impurity entropy 

The impurity entropy can be straightforwardly deter- 
mined from a perturbative expansion of the impurity part 
of the thermodynamic potential. The lowest correction 



to the tree-level value 5*11 



2r In 2 is of order Uq . 



The 



corresponding contribution to the thermodynamic poten- 
tial. Fig. is given by 



dTG}{T)G}{- 



■r). 



(52) 



Here Gf{T) is the fourier-transformed Green's function 
Gf ()26|l in the presence of an UV cutoff A. Taking the 
temperature derivative we can write the entropy result 



A5i 



imp 



2At 



dr G}{t)G}{~ 



(53) 



where Gf — AoGf. The prefactor can be expressed in 
terms of the renormalized coupling and will be pro- 
portional to e at the fixed point. Thus, to leading order 
the square bracket term can be evaluated at r = i, and 
is expected to be a universal, finite number in the limit of 
A — > oo. Unfortunately, we were not able to analytically 
prove the convergence of the integral as A — > oo. We have 
therefore resorted to a numerical evaluation for finite A 



A^i,„p = 0.159y , 



(54) 



Adding the tree-level result and the correction we ob- 
tain 



5imp = ln2 + (0.03±0.005) 



1 



0{e 



3/2 N 



(55) 



The two expansions for the entropy of the SSC fixed 
point, (|22|) and H55|) . predict a small positive correction 
to In 2 for < r < i , as shown in Fig. [TJ consistent 
with the notion that entropy should decrease under RG 
flow^i. Results from NRG (Ref.H as well as ours) show 
that the deviation from In 2 is tiny for all < r < i, 
which is in principle consistent with the analytical re- 
sults. [Although we identified signatures of a correc- 
tion to 5i,np = In 2 in the small-r regime, the accuracy 
of our NRG procedure was insufficient to determine the 
prefactoi»2&.] 



4- Conduction electron T matrix 

In the Anderson model formulation of the impurity 
problem, the conduction electron T matrix is directly 
proportional to the physical / electron propagator. As 
shown in Appendix ^ there are no singular propagator 
renormalizations in the present problem, Z = 1 H40|l. 
Thus, from Eq. H26(l we have the exact result for the ex- 
ponent of the T matrix, 



lmT{Lu) oc \lu\ 



(56) 



which holds at the SCR and SCR' fixed points, and (triv- 
ially) also at the SSC fixed point. A result similar to 
Eq. (|56|l was derived^^ for the SCR fixed point within 
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the small-r expansion of Sec. IIIII in agreement with the 
notion that both expansions describe the same critical 
fixed point. 



V. MAXIMALLY PARTICLE-HOLE 
ASYMMETRIC ANDERSON MODEL 

In this section we consider a different limit of the 
pseudogap Anderson model namely the model with 
maximal p-h asymmetry. This means that one of the 
four possible impurity states will be excluded: with 
C/o — oo (keeping eq finite) this is the doubly occu- 
pied one, whereas with Eq — > — oo (keeping Uq finite) 
the empty state is excluded. As these two situation are 
related by the p-h transformation (0), we will formulate 
the following for [/q — > oo. The Hamiltonian can be writ- 
ten as 



H = eo\<j){<T\+Vo[\(j){e\c^{0) + h.c: 



dk |fc|'' kcl^Cka 



(57) 



where | t), | j), and |e) represent the three allowed im- 
purity states. 

We shall show that this infinite- 1/ pseudogap Anderson 
model has a phase transition between a free-moment and 
a Kondo-screened phase, which is accessible to perturba- 
tive RG techniques near r — 1. In particular, r = 1 plays 
the role of the upper-critical dimension. Furthermore, 
we will argue in Sees. IVII and IVIII that the transitions 
in both the Anderson and the Kondo model with finite 
p-h asymmetry are in the same universality class, and 
are described by the RG presented below, provided that 
r > r*. 



A. Trivial fixed points 

For vanishing hybridization Vb, the maximally asym- 
metric Anderson model features three trivial fixed points: 
for eo < the ground state is the spinful doublet (LM) 
with In 2 entropy. For Eq > we have an empty-state 
singlet, which we can identify with the ASC state of the 
Kondo model (see below). The doubly occupied singlet 
state (labelled ASC') is related to ASC by the p-h trans- 
formation I^J. For Eq = we have three degenerate im- 
purity states, we refer to this as the valence-fluctuation 
fixed point (VFl), with entropy In 3. The impurity spin 
susceptibilities are 



Tx 



imp 



1/4 LM 
1/6 VFl 
ASC 



(58) 



Again, the conduction electron phase shift is zero at these 
fixed points due to the vanishing fixed point value of the 
hybridization. The hybridization term, Vq, is irrelevant 
at the ASC fixed point for all r, and irrelevant at LM for 
r > 0. 



B. Upper-critical dimension: Expansion around 
the valence-fluctuation fixed point 

In the following we perform an expansion around the 
VFl fixed point, i.e., around Eq = 0, Vq — 0. This will 
give access to the properties of the ASC fixed point, i.e., 
a critical fixed point different from the one accessed by 
the RG calculations in Sees. Illll and IIVI 

To represent the three impurity states in the infinite- 1/ 
Anderson model it is useful to introduce auxiliary fields 
for pseudo-particles, &s (for the empty-state singlet) and 
fa (for the spinful doublet). The required Hilbert space 
constraint 6J&s+/1/(t = Q — 1 will be implemented using 
a chemical potential Aq ^ oo, such that observables (O) 
have to be calculated according to^^iffi 



(O) = lim 



(59) 



where (. . .)ao denotes the thermal expectation value cal- 
culated using pseudo-particles in the presence of the 
chemical potential Aq. Clearly, in the limit Aq — > oo 
the term {Q)\o represents the partition function of the 
physical sector of the Hilbert space times exp(— Ao/3). As 
detailed in Ref. ff^, both numerator and denominator of 
Eq. I|59|l have to be expanded in the non-linear couplings 
to the required order when calculating observables; how- 
ever, the denominator does typically not develop loga- 
rithmic singularities at the marginal dimension. 

Furthermore, we need to introduce chemical-potential 
counter-terms which cancels the shift of the critical point 
occurring in perturbation theory upon taking the limit of 
infinite UV cutoff. Technically, this shift arises from the 
real parts of the self-energies of the bs and fa- particles. 
We introduce the counter-terms as additional chemical 
potential for the auxiliary particles. 



SXtblbs, SXfflf 



(60) 



The SXb.f have to be determined order by order in an 
expansion in Vq. Note that counter-term contributions 
in observables in general enter both numerator and de- 
nominator in Eq. H59|l . 

The model (|57() can then be written in the path integral 
form 



dr fa{dr - Sq - Xo - 6Xf)fa 

+ bs{dr - Ao - 6Xb)bs 

+ Vo (fabsCaiO) + C.C.) 



dk \k\'' Ckaidr - k)cka 



(61) 



where Aq is the chemical potential enforcing the con- 
straint exactly, and the rest of the notation is as above. 
The counter-terms ()6Uf) are determined from the real 
parts of the self-energies of both the fa and bg particles 
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at zero temperature according to 

6\f = ReT.f{Xo + eo + SXf,T^O), 
6Xb = ReSb(Ao + (5Afc,T = 0), 



(62) 



note that these real parts diverge linearly with the UV 
cutoff A. 

The model ()61|l shows a transition driven by variation 
of So for finite values of Vq. Tree level scaling analysis 
shows that 



dim[yo] 



1 



r . 



(63) 



This establishes the role of r = 1 as upper-critical dimen- 
sion where Vq is marginal. 

We now proceed with an RG analysis of l|61|l which 
will allow to determine the critical properties for ?' ^ 1 
- a brief account on this appeared in Ref. |l^ The RG 
will treat the auxiliary fields fa and 6s as usual particles 
and consequently determine their propagator renormal- 
izations, anomalous dimensions etc., but all diagrams are 
evaluated with Aq — > oo which ensures the non-trivial 
character of the expansion. Renormalized fields and di- 
mensionless couplings are introduced according to 



fa 
bs 

Vo 



Zf fua 

. V 



(64) 
(65) 

(66) 



No renormalizations are needed for the bulk fermions as 
their self-interaction is assumed to be irrelevant in the 
RG sense. 

In contrast to the field theories analyzed in Sees. IIII Al 
and IIV CI (where the non-linear coupling is used to tune 
the system through the phase transition) , the theory H61|) 
contains two parameters, namely the tuning parameter eg 
and the non- linear coupling Vq. The RG is conveniently 
performed at criticality, i.e., we assume that Eq is tuned 
to the critical line, and RG is done for the coupling v - 
this naturally results in an infrared stable fixed point. To 
two-loop order we obtain the following RG beta function: 



— —rv -\- —V 



(67) 



with the calculation given in Appendix[ni Generally, the 
higher-order corrections to the one-loop result appear to 
be small in the present expansion. One can also consider 
the flow away from criticality, i.e., the flow of the tun- 
ing parameter e, either using S'^ insertions in the field- 
theoretic formulation or explicitly within momentum- 
shell RG. The resulting correlation length exponent is 
in Eq. (IZU below. 

The structure of the above RG is very similar to the 
one of the (4 — e) expansion of the 0'* model for mag- 
nets, where the mass term drives the transition, and the 
non-linear coupling has a non-trivial stable fixed point at 



criticality below four dimensions. Thus, the fixed point 
with finite v (|68|l corresponds to the Wilson-Fisher fixed 
point, whereas w = is the analogue of the Gaussian 
fixed point in the magnetic context, see Fig.|21 The pa- 
rameters V and e play the role of the interaction and the 
mass, respectively. 



C. r* < r < 1 

For r < 1 the trivial fixed point u* = is unstable, 
and the critical properties are instead controlled by an 
interacting fixed point at 



(68) 



At this asymmetric critical fixed point (ACR), we find 
anomalous field dimensions T]b = 2u*^ -I- 2u**, rjf — 
The resulting RG flow diagram is shown 
in Fig. The ACR fixed point (Fig. |2b) shifts to 
larger values of f |£*| with decreasing r, and the ex- 
pansion can be expected to break down for small r. The 
numerical results of Ref. Q; show that this is the case at 
r* « 0.375, where ACR merges with SCR, and p-h sym- 
metry is dynamically restored. Then, Eq. H61I) together 
with an expansion in v, f yields the correct description 
of the critical properties for 0.375 « r* < r < 1. Using 
NRG we have numerically confirmed this expectation, 
i.e., the properties of the critical fixed point of the model 
(|61|l vary continuously as function of r for r* < r < 1. 



D. 



r = 1 



For all r > 1, a phase transition still occurs in the 
asymmetric Anderson and Kondo models, but it is con- 
trolled by the non-interacting VFl fixed point at v = e = 
0. For the marginal case r — 1, i.e., at the upper-critical 
dimension, we expect logarithmic fiow. In the following 
we explicitly keep the UV cutoff A, and discuss RG un- 
der cutoff reduction A AA. We restrict ourselves to 
criticality, where the RG beta function to one-loop order 
is 



dv 
din A 



(69) 



i.e., the hybridization is marginally irrelevant. The RG 
equation can be integrated to give 



1 -3Uo2lnA 



(70) 



with v{X = 1) = Vq. This result will be used below to 
determine logarithmic corrections for a number of ob- 
servables. 
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E. 



r > 1 



For bath exponents r > 1 the couphng Vq in the theory 
(|61|l is irrelevant in the RG sense. The critical system 
flows to the VFl fixed point, Fig. |2|3, and the transition 
becomes a level crossing with perturbative corrections. 

Observables can then be obtained by straightforward 
perturbation theory. Consider, e.g., the boson self-energy 
(Fig. HI; below): 



1 



-k ii^„-iLJ„-Ao-£o ' 



the expression for the ferniion self-energy is similar. At 
the transition, Eq — 0, the self-energies show threshold 
behavior at T = 0, — Iml]y(u; + irj)/n cx V^^(I''"0((D) with 
u) = to — Xq. The low-energy behavior of the fa propaga- 
tor follows as 

-lmGf{Cj + ir])/TT = (1 - A)5{u}) + B\ujY-'^<d{oj) (71) 

with A, i? oc . The propagator has a similar form - 
we will use these results below to explicitly calculate the 
local susceptibility. 



a) 



c) 



d) 





e) 



f) 



^ ^ 



FIG. 8: Feynman diagrams for the infinite- (7 Anderson 
model. Full/ wiggly /dashed lines denote fcr/ba/ca propaga- 
tors, the cross is the counter-term 16011 . a) Bare interaction 
vertex Vq. b) Vq contribution to the partition function, which 
also appears in the denominator of Eq. 15911 . c),d) Diagrams 
entering the local susceptibility, Ximp.imp, to order Vq ■ e) Vo 
contribution to Xu.imp- f) contributions to Xu.u- 



F. Observables near criticality 

We start with the correlation length exponent, v, of the 
asymmetric critical fixed point. In the field-theoretic RG 
scheme it has to be determined via composite operator 
insertions into the action Ht)l|) . which take the role of 
mass terms driving the system away from criticality. The 
lowest-order result for v is 



1. Local susceptibility 

The anomalous exponent i]^ associated with the local 
susceptibility is calculated as above by determining the 
Xioc renormalization factor, Z^, using minimal subtrac- 
tion of poles, and then employing 



Vx 



av 



(73) 



i = r + 0{f^) (r < 1) , 

V 



(72) 



with details of the derivation given in Appendix El For 
r > 1 the transition is a level crossing, formally v — \. 

In the calculation of observables like susceptibilities 
etc. the UV behavior depends on whether we are above 
or below the upper-critical dimension r — \. For r > 1 
the cutoff A has to be kept explicitly, as integrals will 
be UV divergent. This also implies that hyperscaling is 
violated, and no LojT scaling in dynamics occurs, as usual 
for a theory above the upper-critical dimension. 

For r < 1 the UV cutoff A can be sent to infinity af- 
ter taking into account the contributions of the counter- 
terms (|60|) . as the remaining integrals are UV conver- 
gent. This is in contrast to the expansions for r ^ 
f Sec, mill and for r < 1/2 (Sec. lIV|l . which are effectively 
both above a Zower-critical dimension, and where inter- 
mediate quantities can display UV divergencies. Notably, 
for all r < 1 expansions the low-energy observables are 
bdly universal, i.e., cutoff-independent, and hyperscaling 
is fulfilled. 



Here we have xioc = l/"-^ at tree level, in contrast to the 
p-h symmetric problem of Sec. IIVI 

The leading contributions to arise from the dia- 
grams in Fig. 13;, see Appendix lUl for details. From the 



expression 



Z^ = l 



we find T]^ in an expansion in r. 



3(1-0 



(r < 1) 



(74) 



(75) 



This result is again in good agreement with numerical 
findinggi, with a comparison given in Fig. |S1 

We now describe the evaluation of the logarithmic cor- 
rections present at r = 1. Using the lowest-order results 
/3(u) = and ri^{v) = 2v^ we can integrate Eq. l(7^ to 
yield 



(76) 
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The explicit scaling relation^-- for xioc at T = reads 



2. Impurity susceptibility 



Xioc{iio,Vo,A) - ^x(AA)-^xioc (^^,i;(A),A = 1^ . 

(77) 

To analyze the frequency dependence of xioc we choose 
A = A* = cj/A and employ together with to 
express the renormalization factor in terms of A. This 
leads to 



xioc(«w,yo, A) 



UJ 



■Xioc(«,'y(A),l). 



To obtain the leading (i.e. multiplicative) logarithms, the 
last term can be approximated by its fixed point value, 
which simply gives a constant - this neglects sub-leading 
additive logarithmic corrections. We finally obtain the 
result 



1 



w|lnu;|2/3 



(r = l) 



(78) 



valid at criticality for uj <^ A. Note that the structure of 
the logarithms in our problem is different from e.g. that 
of the Hertz-Millis theory at the upper-critical dimension 
studied in Ref. .,32- In our problem both v and flow to 
zero according to Eq. (|76|l : whereas in Ref. l3^ the renor- 
malization factor Z tends to a non-universal constant as 
A — > - this leads to the absence of multiplicative loga- 
rithms. 

Above the upper-critical dimension, r > 1, we have 
simply rix = and thus xioc oc 1/T or c>c l/iv. We 
shall explicitly demonstrate the calculation of the local 
impurity susceptibility using bare perturbation theory. 
To lowest non-trivial order, xioc is given by the convo- 
lution of two fa propagators (|71|l . calculated with the 
self-energy to second order in Vq. Note that no vertex 
corrections occur to this order due to the structure of 
the interaction. When calculating xioc the T — > limit 
has to be taken with care, as the exponentially small tail 
in ImG/(<D) at < contributes to xioc- One obtains 
for the imaginary part Xioc ^'^^ T — > and 1 < r < 2: 



The evaluation of the impurity susceptibility requires 
the summation of the diagrams in Fig. |Ht-f, appearing 
in the numerator of the corresponding Eq. (|59|l . and also 
a careful treatment of the denominator, as we are inter- 
ested in terms being non-singular as function of r. 

The diagrams are conveniently evaluated in imaginary 
time, e.g., the first correction to Ximp,imp in Fig.|Sl: gives 



^ Jo Jo Jo 



(80) 



where the limit Aq — > oo has been taken in the fa- and 
bg propagators. Gco is the conduction electron Green's 
function at the impurity site, which contains a momen- 
tum integral. The other diagrams in Fig. |Slcan be writ- 
ten down similarly. Performing the r integrals first, one 
obtains 



(8c) 
(8d) 
(8e) 

(8f) 



1 

T 



dk 



fc3 



2k 



4 + — tanh — 
2T 



YlA 

rp2 



dkk 



k^ 



^ 2 [3 -h cosh(fc /T)] fc/T - 4 sinh(fc /T) 



A:3[l + cosh(fc/T)] 



4 tanh — - 
2T 



2k k 1 ^ \ 1 2 ^ 

\ tanh — cosh — 

T 2T 2T 



(81) 



where all terms have to be multiplied by e~^°^ /2. Fig.lSJi 
are the contributions from the counter-terms (|60|) , which 
evaluate to (5Af, = 25A/ = 2V§A. The denominator, 
{Q)\o ' being 3e~^°^ to zeroth order in Vo, receives correc- 
tions from the diagram in Fig. [SJa and from the counter- 
terms (EO)), resulting in 



(Q)ao 



dk tanh — , 

kT 2r ' 



(82) 



1 - 2A S{oj)uj B 



6 



T 



r-2 



to be multiplied with e Now we are in the position 
sgn(cj) (79) to collect all contributions to Ximp to second order in Vq: 



This shows that lu/T scaling is violated, as the finite 
frequency part obeys xioc('^ ^ T) cc lu^"'^, but xioc(w <C 
T) (X T~\ 

Moving away from criticality, one of the fa, bg propa- 
gators looses its d{uj — Xo) contribution - for e > {t > 0) 
this is G/^ , indicating that free- moment behavior is ab- 
sent in this regime. Consequently, the zero-temperature 
static local susceptibility is finite, but diverges upon ap- 
proaching the critical point according to xioc oc t^~^. For 
t < 0, Xioc oc 1/T. Thus, the order parameter mioc jumps 
at the transition point for r > 1. 



Axi> 



6T 



^^"^ dx sinh x 



^ 
^/^dxx^ 



3x 



1 + cosh X 

sinh a; 
1 + cosh ; 



(83) 



- 1 



Note that intermediate terms of the form T^^ have can- 
celled in numerator and denominator of Eq. (|59|l . For 
r < 1 the momentum integrals are UV convergent, and 
do not develop poles in f, i.e., the poles present in the 
Ximp, imp diagrams have been cancelled by contributions 
from Fig. |Sl3/f. In the low-energy limit, the product 
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VgT~^ approaches a universal value. Thus Ximp has in- 
deed Curie form, with a fully universal prefactor depend- 
ing on r only. For r > 1 the integrals require an explicit 
UV cutoff, but no correction to the Curie term arises, as 
(Fq^T-^) vanishes as T ^ 0. 

Performing the integral for r < 1 and expressing the 
result in terms of the renormalized coupling v, the impu- 
rity susceptibility reads 



In 2 
~9~ 



(84) 



With the value of the coupling at the ACR fixed point 
(|68|l we finally find, to leading order in (1 — r), 



i- 0.02988(1 -r) + 0(r2) < i) 
^ ('^ > 1) 



(85) 



to be compared with the numerical results in Fig. El 



3. Impurity entropy 

The impurity contribution to the entropy can be de- 
rived from the free energy as above. At the VFl fixed 
point the entropy is S'imp = In 3, and the lowest-order 
correction is computed by expanding the free energy in 
V. Note that this correction vanishes for r > 1, as w* = 
there. 

The calculation of the impurity entropy in the presence 
of a constraint for pscudoparticles has been discussed in 
Appendix C of Ref. IM The limit Aq — !■ oo suppresses 
all contributions from the unphysical part of the Hilbert 
space, in particular disconnected diagrams in the par- 
tition function. Remarkably, this leads to the appear- 
ance of disconnected diagrams in higher-order terms of 
the expansion for the thermodynamic potential fl. The 
second-order diagram for Q, shown in Fig. evaluates 
to 



As expected, the entropy correction is fully universal and 
finite in the limit T — > 0. Inserting the fixed point value 
of the coupling v into Eq. 1)88(1 . we find the impurity 
entropy as 



5'imp — 



In3-^(l-r)-|-0(f2) (r < 1) 



In 3 



(r>l) 



(89) 



have determined S'imp for r values be- 



NRG calculations^ 

low unity; a comparison is shown in Fig.EI As with most 
of the observables obtained within the (1 — r) expansion 
the agreement of the lowest-order result with numerics is 
surprisingly good even for r values well away from unity, 
indicating that higher-loop corrections are small. 

The results for the RG flow and the impurity entropy 
have an interesting corollary: Stability analysis shows 



that for 



< r < 



the RG flow at criticality and small 



p-h asymmetry is from SCR to ACR. The entropy of SCR 
approaches In 2 as r ^ i . Both NRG and the above ex- 
pansion indicate that the entropy of ACR is larger than 
In 2 for r< i, Fig.[7| Thus we have 5imp,ACR > Simp^CR 
for r* < r < i, i.e., the impurity part of the entropy 
increases under RG fiow (see also Fig. below), in con- 
tradiction to the so-called g-theoremSi! (This is not a 
fundamental problem, as the present model has effec- 
tively long-ranged interactions, and is not conformally 
invariant, such that the proof of the ^-theorem does not 
apply.) 



4- Conduction electron T matrix 

The T matrix in the Anderson model is given by 
T{uj) = VQGf{uj). The physical / propagator G/ is a 
convolution of the auxiliary /„■ and bg propagators, i.e., 
the propagator of the composite operator {f^bg). The 
anomalous exponent is obtained from 



din Zt 



dv 



(90) 



An- 



dk 1^ tanh — . (86) 
k 2T ^ ' 



There is also a contribution to fiimp from the counter- 
terms, but this is temperature-independent and does 
not contribute to the entropy. Taking the temperature 
derivative we find: 



AS; 



imp 



3 



2r 



A/T 



dxlxl^ cosh' 



-A/T 



-2 

2 



(87) 



This integral can be performed in the limit of infinite UV 
cutoff. [In contrast, the integral in H8t)|l is UV divergent, 
as discussed in Sec. ^) Expressing the result in terms 
of the renormalized hybridization and taking the limit 
r — > 1 we have 



In 3 



;ln2 



(88) 



As in Ref. |23 we are able to determine an exact re- 
sult for the anomalous exponent, valid to all orders in 
perturbation theory. The argument is based on the dia- 
grammatic structure of the T matrix, namely the relevant 
diagrams can be completely constructed from full v in- 
teraction vertices and full f/b propagators^^. This leads 
to the relation between Z factors 



ZfZb 



(91) 



This equation can be plugged into (|66|l . Taking the 
logarithmic derivative at fixed bare coupling and using 
/3{v)/v = at any fixed point with finite v*, one obtains 
the exact result 



2f 



ImT(aj) cx 



(92) 



Whereas Eq. H92|l applies to the p-h asymmetric fixed 
point (ACR), and is valid for r* < r < 1, the results of 
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Sec UVHl and Ref.^ have established the same critical 
behavior for the T matrix for the symmetric fixed point 
(SCR) for < r < i. Thus, we conclude that all critical 
fixed points for < r < 1 in the pseudogap Anderson and 
Kondo models display a T matrix behavior of lmT{Lu) oc 

The logarithmic correction to the T matrix at r = 1 
are evaluated in a manner similar to the one for the local 
susceptibility above. With (3{v) = and riT{v) = Sv"^ 
- note that jjt = Vf + Vb only holds at one-loop level 
because Z„ = 1 at this order - we can integrate (Pn|l to 
find 



Vo 



(93) 



With the general scaling relation 



T{tLo, Vo, A) = Zt (AA)-ir ( — , v{X), A = 1 ) (94) 



and lf7n|l we find 



T{i,v{X),l). 



As above, the leading logarithms are of multiplicative 
character, with the final result T{uj) oc l/(w| lnw|). This 
implies for the T matrix spectral density the behavior 



Imr(ci;) cx 



1 



w I In cj I ^ 



(r = l). 



(95) 



Above the upper-critical dimension, r > 1, we have 
again r/T = and ImT(ci;) cx S{llj). 



VI. GENERAL PARTICLE-HOLE ASYMMETRY 

Here we comment on the general case of finite p-h 
asymmetry. Starting with the trivial fixed points, tree- 
level power counting shows that LM is stable w.r.t. p-h 
asymmetry, with a scaling dimension of — r. In contrast, 
SSC is unstable towards ASC, and p-h asymmetry grows 
near SSC with a scaling dimension of r. Finally, at the 
free-impurity fixed point (FImp) p-h asymmetry grows 
under RG with a scaling dimension of unity and the sys- 
tems flow towards VFl. VFl itself is stable w.r.t. a devi- 
ation from maximal p-h asymmetry. 

In the r ranges where the RG expansions of this paper 
are perturbatively controlled, we can immediately come 
to conclusions about the stability of the critical fixed 
points: SCR will be stable w.r.t. finite p-h asymmetry 
for r close to 0, but unstable for r close to i. Similarly, 
for r close to or larger than unity, ACR is stable w.r.t. a 
deviation from maximal p-h asymmetry, in other words it 
is safe to discard one of the two impurity states, |d) or |e), 
for the discussion of the critical properties, and to work 
with the perturbative expansion of Sec. One has to 
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jQ-32 jQ-24 jQ-20 jq.16 ^q-12 jq-4 jqO 

Temperature T 

FIG. 9: NRG results for a slightly p-h asymmetric Anderson 
model at r = 0.45 close to criticality (ACR). Parameter values 
are a hybridization strength ttVop{0) — 1, £o = —0.5, and 
Uo = 0.9949, 0.99493 (dashed), 0.9949213, 0.9949215 (solid). 
The thin lines are for a p-h symmetric model very close to 
criticality (SCR), with eo = -0.487, Uo = 0.974. The NRG 
runs were done using a discretization parameter of A = 9, 
keeping TVs ~ 650 levels. The two-stage flow described in the 
text can be clearly seen, i.e., there is a small energy scale 
T* where the system flows from ACR to ASC or LM (solid: 



10" 



dashed: T* 



10" 



and a larger scale Tacr 



T* 

10~^^ where the system flows from SCR to ACR. a) Impurity 
susceptibility Tximp . b) Impurity entropy Simp. The behavior 
of the entropy illustrates the point made in Sec IV F 51 Simp(r) 
decreases as function of T, i.e., increases along the RG flow, 
for T* <T < Tacr due to S'lmp.ACR > Simp.scR. This "uphiU 
flow" does not violate thermodynamic stability criteria, as the 
total entropy (impurity plus bath) of the system still decreases 
under RG. 



keep in mind that ACR moves towards smaller effective 
p-h asymmetry (larger values of |e*|) upon decreasing r, 
and the perturbative expansion around VFl breaks down 
as r ^ r*^ where p-h symmetry is dynamically restored. 

The numerics of Ref. gives no indications for addi- 
tional fixed points in the case of finite p-h asymmetry 
which would not be present in the maximally asymmet- 
ric model; our RG results are consistent with this. Thus, 
the critical properties of a pseudogap Kondo or Anderson 
model with general p-h asymmetry are always identical 
to the ones of the maximally asymmetric model of Sec.lVl 

For the RG treatment of a model with general p- 
h asymmetry close to criticality one has to envision a 
two-step RG procedure, as usual in problems with dif- 
ferent energy scales. Suppose we start near criticality 
from small hybridization and small p-h asymmetry. For 
r > 7'* , the scaling dimension of p-h symmetry breaking 
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term is largest: at tree level we have dini[2eo + Uq] = 1 
and dim[Vo] = (1 — r)/2. Thus the initial model param- 
eters flow towards large p-h asymmetry first. This flow 
drives the system into the regime described by the maxi- 
mally p-h asymmetric Anderson model of Sec.0 thereby 
renormalizing the parameters e and v. Then, the RG of 
Sec. Intakes over, with a flow towards ACR, and finally to 
one of the two stable phases (if the system is not exactly 
at criticality) . Thus, a low-energy scale T* and a higher 
scale Tacr exist, where Tacr characterizes the approach 
of the effective p-h asymmetry towards the ACR fixed 
point. Note that for r > 2 the fastest flow is the one 
of V to zero, i.e., the system quickly approaches the VFl 
fixed point. The described behavior is nicely borne out 
by NRG calculations for the Anderson model, see Fig. IHl 

Finally, for r < r* the initial fiow is dominated by the 
decrease of p-h asymmetry - this cannot be captured by 
our Anderson model RG, but is contained in the Kondo 
treatment ~ and the system approaches SCR, before it 
finally departs to one of the two stable phases. If the 
system is on the strong-coupling side of the transition 
(t > 0), then the behavior near SCR is multicritical (see 
also Fig. 16a of Ref.yi): two low-energy scales exist which 
describe the departure of the fiow from SCR, namely T* 
(for the deviation of j from j*) and a lower scale Tasc 
(for the subsequent growth of p-h asymmetry when flow- 
ing towards the ASC flxed point). 

We note that a recent investigation of the pseudogap 
Anderson model^'^ using the local- moment approach^ has 
found indications of a line of critical flxed points in the p- 
h asymmetric case, parametrized by p-h asymmetry. We 
believe that this is an artifact of the employed approxi- 
mation scheme, as (i) NRG calculations strongly hint to- 
wards a single asymmetric critical fixed point (for flxed 
r), i.e., the fixed-point level spectrum at criticality does 
not depend on the initial p-h asymmetry (e.g. the ratio 
Uo/eq), and (ii) our analytical RG shows that, at least 
near r = 1, the scaling dimension of p-h asymmetry is 
largest, and the sketched two-step RG (which directly 
leads to consider an infinite-f/ Anderson model to de- 
scribe the critical behavior) is a controlled approach for 
arbitrary initial p-h asymmetry. 



VII. RELATION BETWEEN ANDERSON AND 
KONDO MODELS 

This section shall highlight the relation between the 
pseudogap Anderson and Kondo impurity models. 

On the one hand, it is well-known that the Ander- 
son model reduces to the Kondo model in the so-called 
Kondo limit, see Sec. II Al This mapping covers the far 
left-hand side of the flow diagrams in Figs. ^ and and 
suggests that the phase transition at small r in the An- 
derson model in this Kondo limit is described by the 
Kondo RG of Sec. Iml 

On the other hand, we have argued that the flow of 
the Kondo model can be naturally understood in terms 



of the variables of the Anderson model. In particular, the 
RG expansions of Secs. lIVl andlVI describing transitions 
of the Anderson model near r — ^ and r = 1, also ap- 
ply to the Kondo model. Clearly, in thinking about this 
"mapping" of the Kondo to the Anderson model one has 
to view the Anderson model as an effective low- energy 
theory of the Kondo model. The impurity states of this 
effective Anderson model are thus many-body states ob- 
tained after integrating out high-energy degrees of free- 
dom from the Kondo model, i.e., dressed impurity states. 

A formal way to obtain an Anderson model from a 
Kondo model is a strong-coupling expansion: The | t) 
and I I) states of the Anderson model are bare impurity 
states, whereas |e) and |d) represent the impurity with 
either a hole or an electron of opposite spin tightly bound 
to the impurity. 

The above can also be made plausible in a slightly dif- 
ferent way, in the following for a p-h asymmetric situa- 
tion: Let us consider the Kondo problem with a hard-gap 
DOS^-: in the presence of p-h asymmetry it shows a flrst- 
order quantum transition, i.e., a level crossing, between a 
Kondo-screened singlet and a spin-i doublet state. As we 
can understand the asymmetric pseudogap DOS as con- 
sisting of an asymmetric high-energy part and a (asymp- 
totically) symmetric low-energy part, we can obtain an 
effective theory for the pseudogap Kondo model by cou- 
pling the above mentioned three (many-body) impurity 
states, obtained by integrating out high-energy degrees 
of freedom from the band, to the remaining low-energy 
part of the conduction electron spectrum. A crucial in- 
gredient is now the p-h asymmetry of the original model. 
It is clear that upon integrating out the high-energy part 
of the bath two many-body singlet states arise, namely 
|e) and |d) as discussed above. Due to the p-h asymme- 
try of the underlying model these two singlet states will 
have very different energies, such that we can discard the 
high-energy state in the low-energy theory. With this 
we directly arrive at an infinite-fj pseudogap Anderson 
model. 

We conclude that the phase transitions of the pseudo- 
gap Anderson and Kondo models are in the same uni- 
versality classes - this is supported by the numerical cal- 
culations of Ref. '6^. For small r the phase transition is 
naturally described in the Kondo language of Sec. lIIII im- 
plying that the critical fixed point of the Anderson model 
is located in the Kondo limit. In contrast, for larger r the 
formulation in terms of the Anderson model provides the 
relevant degrees of freedom, in other words, both spinful 
and spinless (many-body) impurity states play a role in 
the critical dynamics. 



VIII. FINITE MAGNETIC FIELD 

Interesting physics obtains in the pseudogap Kondo 
and Anderson models in the presence of a finite magnetic 
field. We concentrate here on the effect of a local field, 
applied to the impurity only. Note that a finite field 
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FIG. 10: Schematic T = phase diagram of the pseudogap 
Kondo model in a local magnetic field B. At B = 0, there 
is a phase transition at Jk = Jc between LM and a strong- 
coupling phase, SSC or ASC. Upon application of a field, an 
unscreened spin becomes strongly polarized (POL) , whereas a 
screened spin in a strong-coupling phase is only weakly polar- 
ized. Increasing the field at Jk > Jc drives a phase transition 
- such a transition is not present in the metallic case, r = 0. 
The phase diagram for a pseudogap Anderson model is simi- 
lar. As discussed in the text, the zero-field critical fixed point 
is always unstable w.r.t. finite B in the RG sense; thus the 
two transitions are in general in difi'erent universality classes 
(which is already clear from symmetry considerations). 



applied to the bulk can modify the low-energy behavior of 
the bath DOS due to Zeeman splitting; then, pseudogap 
Kondo physics survives only for energy scales above the 
Zeeman energy. 

For T = and a metallic density of states, r = 0, a 
local field B only leads to a crossover as function oiB/TK, 
and all observables evolve smoothly from B = (screened 
spin) to large B (polarized spin). In contrast, for any 
r > and t > (i.e. Jk > Jc in the Kondo model) there 
will be a zero-temperature phase transition as function of 
the local field between one phase with weak and one with 
strong impurity spin polarization, see Fig. 1101 We briefly 
discuss this phase transition in the following - it turns 
out that the variables of the Anderson model provide a 
natural language to analyze the problem. Remarkably, 
all ingredients needed for the discussion of the critical 
properties have already been calculated in the previous 
sections. 

Strategically, we first discuss the decoupled impurity, 
and then analyze the modifications arising from the pres- 
ence of the hybridization term. For the decoupled impu- 
rity, the local field is trivially a relevant perturbation, 
with scaling dimension unity. Thus, in the low-energy 
limit the field is effectively infinite, and one of the two 
impurity states | t): I i) can be discarded (we will discard 
I I) in what follows). 



A. Asymmetric Anderson model 

In the p-h asymmetric case two impurity states have 
to be considered, namely | t) and |e). We discuss a level 



crossing transition of these two, tunable e.g. by vary- 
ing the magnetic field, and being coupled to conduction 
electrons. At the transition we are left with a (spinless) 
resonant level model, and the analysis for a pseudogap 
host density of states is in Sec. IIVBI In particular, for 
r > 1 the phase transition is a level crossing with per- 
turbative corrections, and for all < r < 1 we have a 
continuous transition with a critical fixed point identical 
to the intermediate-coupling fixed point of the resonant 
level model, Eq. ||2U. Thus, the properties of the tran- 
sition evolve smoothly as function of r for < r < 1, in 
contrast to the zero- field situation (!). 



B. Symmetric Anderson model 

In the presence of p-h symmetry and magnetic field, 
the decoupled impurity has three low-energy states: | t), 
|e), |d). The resulting level-crossing transition is tech- 
nically identical to the one of the zero-field infinite-t/ 
Anderson model of Sec. El the two situations can be 
mapped onto each other via the p-h transformation iPJ. 
From Sec. we can read off the properties of the field- 
tuned phase transition: For r > 1 we have again a level 
crossing with perturbative corrections. A RG expansion 
around the level-crossing fixed point can be used to cal- 
culate the critical properties below r = 1. This expansion 
describes the physics for all r > r* . 

It is interesting to ask what happens for r < r* . At r* 
the broken level symmetry is dynamically restored at the 
critical fixed point, meaning that for r < r* the tran- 
sition is controlled by a zero- field critical fixed point. 
From Sees. Illll and IIVI we know that for r < ^ such an 
Anderson model transition with four degenerate impu- 
rity states is the same as described by the weak-coupling 
RG in the Kondo language. However, the mapping Q 
shows that we have to consider a Kondo model of a 
charge pseudospin here, i.e., the corresponding Schrieffer- 
Wolff transformation will project out the | t) and | |) 
states. In other words, in a p-h symmetric situation with 
< r < r* , where the zero-field transition is controlled 
by the SCR fixed point, the finite-field transition will be 
asymptotically controlled by the SCR' fixed point! 



C. Kondo model 

The above statement is consistent with a simple gener- 
alization of the weak-coupling RG of Sec. Illll to the case 
of finite field. One easily finds that the field is a rele- 
vant perturbation at the SCR fixed point with a scaling 
dimension of unity, and the finite-field transition is not 
accessible in a description using the Kondo model spin 
variables. (In contrast, the magnetic field is irrelevant at 
SCR'.) 



21 



D. Symmetries and pseudospin field 

In the Anderson model language, the local field lifts the 
degeneracy of the magnetic doublet | I i)- It breaks 
the SU(2) symmetry in the spin sector, and the effect is 
clearly independent of the field direction. 

As discussed in Sec. II Al the p-h symmetric Anderson 
model display SU(2) symmetry also in the charge sec- 
tor (charge pseudospin). This can be broken by choos- 
ing Uq — 2£oi corresponding to a pseudospin field in 
z direction. Clearly, a pseudospin field in perpendicu- 
lar direction will have the same effect - this field corre- 
sponds to a local pairing field coupling to f^fi- Thus, 
the SU(2) pseudospin symmetry shows that the effect 
of s-wave pairing correlations on the Kondo problem is 
similar to a particle- hole asymmetry, as argued in Ref.lssI 
using a completely different approach. 

IX. CONCLUSIONS 

In this paper, we have analyzed the zero-temperature 
phase transitions in the pseudogap Kondo problem, char- 
acterized by a bath density of states p{u!) cx \uj\^. Using 
various perturbative RG expansions, formulated in the 
degrees of freedom of either the Kondo or the Ander- 
son impurity model, we have developed critical theories 
for the phase transitions in both the p-h symmetric and 
asymmetric cases. We have established that r — and 
r = \ pl^y the role of two lower-critical dimensions in the 
p-h symmetric case, i.e., the non-trivial phase transition 
disappears both as r — > 0^ and — > ^ with diverging 
correlation length exponent. In contrast, r = 1 is an 
upper-critical dimension for the p-h asymmetric model. 
The transitions for < r < 1 are described by inter- 
acting field theories with universal local-moment fluctu- 
ations and strong hyperscaling properties including uj/T 
scaling of dynamical quantities. In contrast, for r> 1 we 
find a level crossing with perturbative corrections, and 
hyperscaling is violated. 

In the p-h symmetric case, we found two different ex- 
pansions, one around r = and one around r = i, to 
describe the same critical fixed point. In the presence of 
p-h asymmetry, a different critical fixed point emerges, 
which can be analyzed in an expansion around r — \. 

Apart from the small-r expansion of Sec. IIIII all our 
theories were formulated using the Anderson model lan- 
guage. This shows that the quantum phase transition 
between a screened and an unscreened moment in Kondo- 
type models can be nicely interpreted by saying that the 
system fiuctuates between "possessing a moment" and 
"possessing no moment" - this is precisely what is de- 
scribed by the effective Anderson model at its valence- 
fiuctuation fixed point. 

We have calculated a number of observables, with re- 
sults being in excellent agreement with numerical data. 
In particular, we have found an exact exponent for the 
conduction electron T matrix, valid for all expansions 



used in this paper. We have also discussed the physics of 
the pseudogap Kondo problem in a local magnetic field, 
where we have shown that - in contrast to the metallic 
Kondo effect - a sharp transition occurs as a function of 
the field. 

Applications of our results include impurity moments 
in unconventional superconductors^iifliSiS^ and other 
pseudogap systems, like e.g. graphite. On the theoret- 
ical side, we expect that the analysis of quantum im- 
purity models using Anderson instead of Kondo model 
variables may be useful in a variety of problems. Thus, 
field theories similar to ours can possibly be constructed 
for other impurity quantum transitions, and will also be 
useful for the study of lattice models in dynamical mean- 
field theory^ and its extensions'^^, where local quantum 
criticality can be captured using effective single-impurity 
models. 
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APPENDIX A: THE NON-INTERACTING 
RESONANT LEVEL MODEL 

Here we provide a few details about the non-interacting 
resonant level model with a pseudogap density of states. 
The model is exactly solvable, and all properties can be 
directly evaluated using the dressed / electron propa- 
gator H26() . The local susceptibility has been quoted in 
Sec. UVBl The T matrix is given by T(w) = V;^Gf{w)] a 
brief discussion of spectral properties can also be found 
in the Appendix of Ref. 

We now sketch the calculation of Tximp and S'imp- The 
diagrams contributing to T^imp are in Fig. 1111 note that 
these contain full / propagators, and no higher-order di- 
agrams appear. Power counting shows that terms of l/T 
form do only arise from Xu,u; both Ximp.imp and Xu,imp are 
less singular. The evaluation of the diagrams in Fig. Illb 
gives 

Xu.u(T) ^^TY^r rffclfcr '""'"^ (Al) 

n ^ ' 
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c) 
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FIG. 11: Feynman diagrams for the spin susceptibility of the 
non-interacting resonant level model of Sec. li V lH Dashed/full 
lines denote / dressed fa propagators, respectively; the full 
dots are Vo vertices, and the open circles are sources, a) 
Ximp.imp. b) Xu,imp. c)xu,u- Note that the displayed diagrams 
give the exact results for the susceptibilities, as all self-energies 
are contained in the dressed / propagator. 

All factors of Vq cancel exactly against the Aq, and we 
obtain a quantity with a universal prefactor. The k inte- 
grals can be performed in the limit of infinite UV cutoff, 
followed by the Matsubara summations. The result is of 
Curie form, with 

XuAT)^-^. (A2) 

Note that we have used the dressed / propagators in 
their low-energy form - due to their singular nature their 
high-energy properties are unimportant for the leading 
low-temperature behavior of the susceptibility. 

The impurity entropy S'imp is directly calculated from 
the full / propagator l|2f>|l . Performing the tr In in the 
partition function gives 

n-nup = -2r^ ln[zsgnK)l^^»r] ■ (A3) 

n 

One can easily see that the T = entropy contributions 
of ln(iaj„) and In \iuJn\ are identical, thus we have 

5in,p = 9t ||-2rr^ln(ic^„)e*""0^j . (A4) 

This is nothing but 2r times the entropy of a free spinless 
fermion, resulting in S'imp = 2rln2. 

Interestingly, the above exact universal results can also 
be obtained in the framework of a RG expansion in the 
hybridization, i.e., in a formal expansion around r = 1. 
Here one exploits the intermediate-coupling nature of 
the fixed point, by calculating observables using renor- 
malized perturbation theory and employing the fixed- 
point value H31() of the coupling v. The necessary dia- 
grams. Fig. 1121 now contain bare (instead of dressed) / 
propagators. Technically, higher-order terms appear in 
this expansion, but upon interpreting the lowest-order 
result of renormalized perturbation theory (e.g. re- 
exponentiating logarithms) these terms are completely 
summed up. Thus, within the RG framework the in- 
formation about the non-trivial / scaling dimension is 



FIG. 12: Diagrams occurring in the perturbative RG for the 
non-interacting resonant level model of Sec. IIV Bl Here, full 
lines are bare fc propagators, a) Bare interaction vertex Vq. 
b) fa- self-energy, c) Diagram for the entropy correction, d) 
Diagrams entering the local susceptibility to order v'^. Within 
the renormalized perturbation theory scheme all corrections 
beyond second order in v vanish. 

contained in the coupling value v* instead of in the prop- 
agator G f . We briefly demonstrate the idea by evaluating 
the impurity entropy. As we expand around v = (the 
FImp flxed point) we obtain a perturbative correction to 
the impurity thermodynamic potential which is of order 
Vq (Fis. 112b). Taking the temperature derivative leads 
to 

A5imp = --^T-^' / dx\xY cosh-2 (A5) 

with f = (1 — r)/2 as above. Performing the integral 
in the limits A cxd and r ^ 1 and introducing the 
renormalized hybridization v yields 

5i,„p = (l-2i;2) ln4. (A6) 

With the fixed point value of the hybridization we obtain 
the result S'imp — 2r In 2 as above. 

APPENDIX B: RG FOR THE INTERACTING 
RESONANT-LEVEL MODEL 

In this appendix we present details of the renormaliza- 
tion group treatment for the symmetric Anderson model 
in the vicinity of r = i. It is based on an expansion 
around the SSC fixed point, i.e., around a non-interacting 
resonant-level model, with the interaction strength being 
the expansion parameter. Thus the expansion is done 
around an intermediate-coupling fixed point (!). 

The starting point is the action (|36|l derived in 
Sec. IIV CI for < r < 1. Importantly, the "bare" (i.e. 
C^o = 0) / propagator in behaves as 1/(^0'^'^) at low 
energies, Eq. H26() . as the conduction electrons have al- 
ready been integrated out. The prefactor will be kept 
explicitly in the renormalization steps. We employ the 
field-theoretic RG scheme in order to determine the flow 
of the coupling renormalized coupling m, introduced in 
Eq. (|39|l . with dimensional regularization and minimal 
subtraction of poles. As will be seen below, the lowest 
non-trivial renormalizations arise at two-loop order. 



23 



a) 



c) 



XXX XIX 
XBC XBC ¥^ 



FIG. 13: Diagrams for the RG of the symmetric Anderson 
impurity model, with notation as in Fig. 0| a) Uq self-energy 
contribution, which vanishes, b) One-loop vertex renormal- 
izations, which vanish due to p-h symmetry, c) Two-loop 
vertex renormalizations - these are the only contributions to 
the beta function 1411 1. 



The needed diagrams arising in the RG treatment are 
displayed in Fig. 1131 A number of observations can be 
made: (i) Hartree diagrams vanish due to the overall p-h 
symmetry of the model, (ii) Explicit calculation shows 
that the one- loop vertex renormalization diagrams, i.e., 
the corrections to u in Fig.lOb. do not develop e poles. 
This also follows from the fact that the beta function 
cannot contain even powers in u due to p-h symmetry, 
(iii) Explicit evaluation shows that the / self-energy up 
to two-loop order. Fig. [T^ . does not contribute singular 
propagator renormalizations. In other words, the field 
renormalization factor is 

Z=l (Bl) 

to two-loop order. As we will argue below, this result 
is exact to all orders in perturbation theory. It implies 
that the / field does not acquire an anomalous dimension, 
r?/ = 0. 

The two-loop diagrams for the renormalization of u 
can be divided into two groups of three each, displayed 
in the two lines of Figs. 113b. The diagrams in each group 
can be easily seen to be equal. Explicitly: 



(13cl) = 



duji 



-\ 2 



G f{iuJi)G f{iuji + w) 



(B2) 



where v is the sum of two external frequencies. Using 
the explicit form of G/ (|26|l . the frequency integral can 
be split according to the absolute value in G/ , and then 
performed directly with UV cutoff sent to infinity. The 
result contains e poles and reads 



(13cl) = Ul 



1 



2 In 4 



(B3) 



The second group of diagrams in Fig.ll3h is slightly more 
complicated. All three of them can be brought into the 
form 

dwi du>2 

X Gf{—iuji + ii')G f{ibJ2 — iiui + iuj) 



(13c4) = "^^0 J ■ 



■Gf{iui)Gf{uj2) 



(B4) 



where the prefactor 2 accounts for the different arrange- 
ment of the vertices, and w is a (third) external frequency. 
The llJ2 integral can be performed directly and leads to 
the same expression that occurred in the brackets in (|B2p . 
The remaining uji integral can also be performed after 
splitting the interval into four parts. Straightforward al- 
gebra, using the identity arcsin y^— iarctanh-^/l — l/x ~ 
7r/2, where x = u/v < 1, leads to 



(13c4) = Ul 



2e 



•(B5) 



Upon adding the contributions fromll3bl-ll3b6. only sin- 
gle poles in e remain - as Z = 1, the cancellation of the 
double poles is a consistency check of our calculation. 
Minimal subtraction of poles at external frequencies set 
to /i yields the renormalization factor for the quartic cou- 
pling: 



Zi = 1 



3(7r - 21n4)u2 



(B6) 



The RG beta function can be evaluated by taking the /i 
derivatives of at fixed bare coupling. 



P{u) = fl 



du 
dfi 



3(7r - 2 In 4) 



(B7) 



which is the result in Eq. I|41() . 

We now provide the proof for Z ~ 1, i.e., no singular 
propagator renormalizations occur in the present prob- 
lem. What is needed is p-h symmetry and the form of 
the bare propagator, G/ oc or G/(r) cx t~^^'^ in 

the low-energy limit for r — ^. We argue that no log- 
arithms arise when evaluating the self-energy at r = i. 
By power counting, all / self-energy diagrams in the r 
domain give a result proportional to t~^^^. Let us dis- 
cuss internal t integrals arising in the evaluation (which 
could produce logs). In the first step, those involve four 
Gf propagators due to the quartic interaction - here it 
is important that tadpoles do not contribute due to p-h 
symmetry. Therefore, the integrand will behave as t~'^, 
and logs cannot occur. In further steps, four propagators 
do not necessarily occur, but the power counting argu- 
ment shows that each internal time still has to occur with 
a power of t~^. This shows that the integrands will al- 
ways be more singular than 1/r, and no logarithms arise. 
Clearly, such an argument cannot be constructed for the 
vertex diagrams, as they behave as 1/t at r = ^, and 
Fourier transformation generically yields an e pole. 

Let us turn to the local susceptibility renormalization 
factor Z^. As described in Sec. II V HI one can either 
determine the composite operator renormalization factor 
Z2 and use Z^ = Z|, or evaluate Z^ directly with the 
help of the perturbative correction to xioc in Fig. In 
both cases, the / bubble diagram is needed. To leading 
order in e it evaluates to 



G/(iw„)G/(«w„ + iu) 



(B8) 
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in the zero-temperature limit. For the correlation func- 
tion Ci^'^^ in Fig.Et) we introduce a renormalization fac- 



tor as follows: G 



(2,1) 



= Z-iZ-^G 



(2,1) 

f,R ■ 



Demanding that 



the renormalized ci^'i^ is free of poles at scale fi, and 



using Z = 1 we find: 



a) 



b) 



< >. > 



Zo 



u 
ire 



(B9) 



which gives Z^ in Eq. H45|l . 



APPENDIX C: RG FOR THE INFINITE-C/ 
ANDERSON MODEL 

The infinite-C/ Anderson model can be analyzed in the 
vicinity of r = 1 using an expansion in the hybridization 
strength. Here we describe details of the field-theoretic 
RG, with the results appearing in Sec.0 Starting point 
is the action H61(l . written in terms of auxiliary fields fa 
and bg. The RG proceeds by perturbatively calculating 
renormalizations to the propagators and vertices appear- 
ing in l|61() . which yields expression for the renormaliza- 
tion factors defined in Eq. Hfififl . 

The relevant diagrams are displayed in Fig. E| Note 
that there is no one-loop contribution to the vertex renor- 
malization. At one-loop order we have the graphs for the 
/ and b self-energies. Evaluation of the first diagram for 
the / self-energy in Fig. 114b . gives: 



dk\k\ 



1 



-k iv — iijJn — Xa 



(CI) 



We work at criticality, i.e., Eq — 0. First the frequency 
sum can be performed, here the Aq — > oo limit is impor- 
tant: it suppresses contributions from positive k. The 
remaining k integral is UV convergent, and we finally 
obtain 
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FIG. 14: Diagrams for the RG of the infinite-[/ Ander- 
son model. Notations are as in Fig. |H| a),b) Self-energies 
for /ct and bs to two-loop order. The crosses denotes the 
counter-terms which cancel the real parts of the self-energy 
insertions, c) Two-loop vertex renormalization; there is no 
contribution to one-loop order, d) Diagrams with mass in- 
sertions (squares), needed to determine the correlation length 
exponent. 



Collecting all expressions yields the renormalization 
factors, defined in Eq. ^<6<6\ . to two- loop order as 
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With these results, we can take the logarithmic fi deriva- 
tive of Eq. at fixed bare coupling, we obtain the beta 
function 



(C5) 



Vo 



as quoted in Eq. H67|) . The anomalous field dimensions 
are obtained from 



In the second line we have expressed the result in 
terms of renormalized quantities, with = ^"^"^{iv — 
Xo)^'^^ Z^ / {Z f Zb) . Demanding cancellation of poles in 
the expressions for the renormalized / Green's function 
at external frequency iv — Xq = fi we obtain the one-loop 
result for the / renormalization factor as 
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The other diagrams in Fig. I14b.b.c are evaluated in a 
similar manner. In the two-loop self-energies, the real 
part of the inner self-energy insertions, which diverges 
linearly with the UV cutoff, is exactly cancelled by the 
corresponding counter-term. 



dlnZ 



Vf 
Vb 



f 



P(^^)'^h^^y2^2v^ (C6) 
dv 



dfi 



dii 



dv 



To determine the flow away from criticality and the 
correlation length exponent, we follow the standard 
scheme^-^ via insertions of composite operators, repre- 
senting mass terms. Physically, only the difference be- 
tween the masses of the / and b auxiliary fields is rel- 
evant, and introducing one type of composite operator 
is sufficient. In the following we work with // inser- 
tions, which acquire a corresponding renormalization fac- 
tor Z2f. To determine Z2/ we consider a correlation 
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function GiJ'i' = ((fot (r)6(r'); (/t/)(r"))). The 
malization mctor Z2 / is then defined through 



APPENDIX D: SPIN ANISOTROPIES 
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Evaluating the diagram in Fig. I14t i and demanding can- 
(2 1) 

cellation of poles in G),y^ gives, to one- loop accuracy, 

Z2/-l + ^. (C9) 
r 

The correlation length exponent can be related to the Z 
factors via the RG equation^^ for the renormalized G^^j,'^^ , 
according to 



1 , ^2/ 

- - 1 = ^— In— ^ 
ly ail Zf 



(CIO) 



Thus far our discussion has been restricted to impurity 
models with full SU(2) spin symmetry. In this appendix 
we show that spin anisotropics, e.g. in the Kondo inter- 
action, are irrelevant at all critical points considered. 

For the weak-coupling RG of Sec. IIII Al formulated in 
the variables of the Kondo model, we introduce Kondo 
couplings Jj_ and Jz for the transverse and longitudinal 
part of the Kondo interaction. The RG equation H18|) , to 
one-loop order, generalizes to 



j1 



(Dl) 



which yields 



1 -,0*2 
— = 1 — 3v 

V 



(Gil) 



The result for v is in Eq. (I72[l . and is of course identical to 
the one obtained in Rcf. 12 using the familiar momentum- 
shell method. 

The local susceptibility exponent is determined via a 
renormalization factor, Z^, for the two-point correlations 
of the impurity spin. The leading diagrams are in Fig.|SJ:; 
the diagrams are identical, with a leading singular 
contribution of 



TtK,' / dT2 / dTlG,o(T2 - Tl) = ^T^' 

2 Jo Jo 4r 

each. Note that the counter-term contributions can be 
ignored here, as they do not develop poles in f . Demand- 
ing cancellation of poles yields 
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Taking the logarithmic derivative w.r.t. /i one obtains 
the 77^ value quoted in Sec. IV Fl 

The renormalization factor for the T matrix (diagrams 
not shown) is obtained as 
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2 f 



(G13) 



to two-loop order. The resulting anomalous exponent tjt 
fulfills the exact equation rjT = 2f . 



For the metallic case, r = 0, there is a line of fixed points 
at j_L = 0, jz<0, representing an unscreened moment. 
In contrast, for r > only SU(2) symmetric fixed points 
survive, namely LM with = J2 = ^^nd SGR with 
j± ~ Jz = f- Thus, the symmetric critical fixed point of 
the Anderson and Kondo models is stable w.r.t. SU(2) 
symmetry breaking. 

Turning to the p-h asymmetric model, we start from 
the infinite-?/ Anderson model (|57|l and introduce a spin 
dependence in the hybrization, Vb| 7^ Vo|. As is easily 
seen via Schrieffer- Wolff transformation, for a generic p- 
h asymmetric conduction band this is equivalent to an 
anisotropic exchange interaction plus a local magnetic 
field. Without loss of generality we restrict the following 
analysis to a p-h symmetric conduction band, where the 
field term is absent. The RG equation (|67|l . to one-loop 
order, generalizes to 



/9(«t) = -^«T + y(2i'T+«i)' 



(D2) 



Apart from the SU(2) symmetric fixed point, u^* = |f, 
there are two other fixed points with — r, vi ~ and 
vj* = f, W| = - these are, however, infrared unstable 
w.r.t. finite W| (w|). Thus, the only stable critical fixed 
point is the SU(2) symmetric one, which corresponds to 
the AGR fixed point analyzed in Sec. 
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